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Knowledge  of  the  energy  balance  and  plasma  motions  in  the 
atmosphere  is  essential  for  the  accurate  prediction  of  atmo¬ 
spheric  effects  on  radio  communications  under  both  normal 
and  disturbed  conditions.  A  substantial  portion  of  the 
variable  content  of  atmospheric  energy  is  stored  in  charged 
particles.  Previous  investigations  have  demonstrated  that 
the  concentrations  of  atomic  ions  are  controlled  by  transport. 
Electric  fields  play  an  important  part  in  determining  transport 
properties.  Section  2  reviews  some  aspects  of  the  photo¬ 
chemistry  of  atmospheric  ions  and  assesses  the  magnitude  of 
transport  effects  on  atomic  ions. 

The  effort  undertaken  under  the  present  investigation  has 
focused  on  formulating  procedures  to  utilize  satellite  measure¬ 
ments  in  order  to  obtain  a  self-consistent  model  of  charged 
particle  motions.  Such  a  formulation  is  presented  in  Section  3. 
Section  4  deals  with  the  simpler  situation  of  one  ionic  species 
in  addition  to  electrons  for  the  purpose  of  illustrating 
several  features  of  the  general  theory.  Results  utilizing 
satellite  data  are  presented  in  Section  5.  It  is  shown, 
generally,  that  the  major  ion  moves  in  response  to  the  sum  of 
the  forces  on  itself  and  on  electrons,  and  that  minor  ions 
follow  the  motion  of  the  major  ion.  It  is  further  shown  that 
the  observed  vertical  plasma  drift  velocity  can  be  accounted 
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for  by  an  eastward  electric  field  of  between  3  and  5  millivolts 
per  meter.  Calculations  based  on  the  assumption  that  horizontal 
variations  are  negligible  cannot  account  for  the  transport 
term  in  the  equation  of  continuity,  probably  because  the  data 
pertains  to  the  (magnetically)  equatorial  region.  Application 
of  the  present  approach  to  midlatitude  data  is  suggested. 

Appendix  A  lists  the  data  used  in  the  calculations. 

Appendix  B  reviews  the  definition  of  momentum  transfer  colli¬ 
sion  frequencies,  and  presents  collision  frequencies  for 
electrons  and  0+  ions.  An  updated  version  of  atomic  and 
molecular  cross  sections  is  presented  in  Appendix  F.  A  pro¬ 
cedure  for  the  calculation  of  the  spectrum  of  secondary  elec¬ 
trons  produced  by  energetic  precipitating  particles  is  given 
in  Section  6,  and  Appendices  C,  D,  and  E. 
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2. 


ION  CONCENTRATIONS  IN  THE  DAYTIME  IONOSPHERE 


The  altitude  profiles  of  ionic  species  in  the  upper 
atmosphere  are  determined  by  photochemical  processes  and  by 
transport.  In  order  to  assess  the  effects  of  transport,  it  is 
necessary  to  have  knowledge  of  the  relative  importance  of  the 
photochemical  effects.  Using  data  from  an  Atmosphere  Explorer 
orbit  (see  Appendix  A  for  details)  we  show  that  the  altitude 
profiles  of  molecular  ions  are  determined  by  photochemical 
processes  under  the  assumption  of  local  equilibrium,  i.e.  that 
transport  effects  are  negligible  for  molecular  ions.  For 
atomic  ions,  on  the  other  hand,  transport  processes  play  the 
dominant  part  in  determining  altitude  profiles,  particularly 
at  the  upper  altitudes. 

Fig.  1  is  a  plot  of  the  concentration  of  the  molecular 
ions  (0+,  N0+,  N*)  versus  altitude.  The  observed  concentrations 
are  in  harmony  with  theoretical  values  calculated  under  the 
assumption  of  local  photochemical  equilibrium.  The  agreement 
between  theory  and  observations  is  due  to  the  short  chemical 
lifetimes  of  the  molecular  ions.  Fig.  2  is  a  plot  of  the 
chemical  lifetime  of  each  ion  vs.  altitude.  The  primary  sink 
for  molecular  ions  is  dissociative  recombination.  Accordingly 
chemical  lifetimes  for  these  ions  tend  to  decrease  with 
increasing  altitude,  following  the  increase  with  altitude  of 
the  electron  concentration.  Under  daytime  conditions,  the 
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sources  are  balanced  by  the  sink  from  dissociative  recombination 
in  the  equation  of  continuity.  At  night,  both  of  these  terms 
are  reduced  significantly,  and  transport  terms  become  of 
comparable  magnitude. 

Fig.  2  further  shows  that  chemical  lifetimes  increase 
rapidly  with  altitude  for  the  atomic  ions.  The  increase  is 
particularly  rapid  for  He+,  for  which  the  only  significant  sink 
is  the  reaction  with  molecular  nitrogen.  Thus  the  chemical 
lifetime  of  He+  is  inversely  proportional  to  the  ^  concentra¬ 
tion,  and  therefore  transport  effects  dominate  the  He+  profile 
at  the  upper  altitudes.  This  is  apparent  in  Fig.  3,  which  shows 
the  observed  and  calculated  altitude  profiles  for  the  atomic 
ions . 

Fig.  3  contains  three  profiles  for  each  of  the  other 
atomic  ions  (0+,  H+,  and  N+) .  At  the  upper  altitudes,  the 
significant  sources  for  these  ions  are  photoionization  and 
charge  transfer  from  other  atomic  ions.  The  solid  line  profile 
for  each  ion  is  a  purely  theoretical  calculation  under  the 
assumption  of  local  photochemical  equilibrium  (i.e.  no  transport 
effects  are  included)  and  the  concentrations  of  the  ions  are 
determined  by  solving  the  coupled  equations  of  continuity. 

Fig.  3  shows  that  these  purely  theoretical  profiles  are 
increasingly  at  variance  with  the  observed  concentrations 
above  250  km,  reflecting  the  neglect  of  transport  on  all  the 
coupled  species.  Fig.  3  also  displays  a  second  theoretical 


profile,  calculated  under  the  assumption  of  photochemical 
equilibrium  but  utilizing  observed  concentrations  for  all 
ions  except  the  ion  of  interest.  In  this  instance  the  agree¬ 
ment  between  calculated  and  observed  values  is  very  good. 

This  latter  result  indicates  that  the  model  of  photochemical 
reactions  utilized  is  an  accurate  one.  It  does  not  follow, 
however,  that  the  atomic  ions  are  in  local  photochemical 
equilibrium.  In  what  follows  we  indicate  why  the  agreement 
between  observed  ion  concentrations  and  those  calculated  by 
the  photochemical  model  (without  transport)  improves  if, 
in  computing  the  concentration  of  a  particular  ion,  the  observed 
concentrations  of  all  other  ions  are  used.  The  improvement  is 
very  substantial  in  the  case  of  0  ,  H  and  N  .  The  improvement 
results  from  the  fact  that  the  observed  concentrations  impli¬ 
citly  contain  the  effects  of  transport.  We  demonstrate  this 
idea  by  considering  two  ions  whose  photochemistry  is  coupled. 

The  equation  of  continuity  states  that  the  rate  of  pro¬ 
duction  minus  the  rate  of  loss  equals  the  divergence  of  the 
flux: 

P  -  L  =  V  •  F  . 

For  the  example  under  consideration,  we  have 


P1  =  Q1 

+  k2X2 

=  («,1+k1)X1 

H  t^ 

L2  =  (i2+k2)X2 

-+ 

P2  =  Q2 

+  klXl 

V  •  F 

2 

2  t2X2 
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where  the  Q's  represent  production  rates  not  involving  the 
ions  of  interest,  the  2,’s  represent  loss  rates  not  involving 
the  ions  of  interest,  the  X's  represent  the  ion  concentrations, 
and  the  t's  are  defined  by  the  equations.  The  equations  of 
continuity  may  then  be  written  as 

Q1  +  k2X2  "  <ei+kl>Xl  =  fcixi 
q2  +  -  (d2+k2)x2  =  t2X2  . 

The  solution  is 

X1  =  f(e2+k2+t2)  Q1  +  k2°25  /  n^2+k2t2) (£l+kl+tl)  "  klk2] 

with  an  analogous  expression  for  X2.  This  is  the  exact  solu¬ 
tion  and  therefore  represents  the  observed  value  of  X^,  If  we 
neglect  transport  (t^  =  t2  =  0)  the  solution  is 

Yx  =  [(<>2+k2)  Ql  +  k2Q2]  /  [  (Jl2+k2)  (Jl1+k1)  -  kxk2] 

where  a  different  symbol  (Y)  has  been  used  to  indicate  that 
this  is  an  approximate  value.  We  compare  this  with  the  exact 
(observed)  value  X^ : 

Yf  =  Xi  {  i  +  [(^2+k2)t1  +  k2t2(x2/x1)]  /  U1£2+£ik2+e2kl)  } 

We  note  first  that  the  difference  is  of  order  t /%.  This  is  a 
ratio  that  increases  rapidly  with  altitude,  since  transport  ef¬ 
fects  increase  with  altitude  whereas  £  decreases  in  proportion 
to  the  neutral  density.  Secondly,  the  difference  involves, 
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in  general,  the  transport  of  both  ions.  Thirdly,  the  differ¬ 
ence  can  be  of  either  sign,  since  the  t's  can  be  of  either 


sign.  Finally,  if  the  ratio  of  concentrations  (X^/X^)  is  very 
large,  the  inaccuracy  in  the  concentration  of  the  minor  ion 
depends  on  the  transport  characteristics  of  the  major  ion. 

This  observation  indicates  that  the  use  of  observed  concentra¬ 
tions  is  particularly  beneficial  to  the  calculation  of  minor 
ion  concentrations. 

Using  observed  concentrations  in  a  pure  photochemical 
calculation  (without  transport)  yields  the  solutions 

Zx  =  Xj_  [1  +  t1/(£1+k1)  ] 

We  have  used  a  different  symbol  (Z)  for  the  computation  that 

neglects  transport  but  uses  observed  concentrations  for  ions 

other  them  the  ion  of  interest.  The  error,  in  this  instance, 

involves  only  the  transport  properties  of  the  ion  of  interest. 

Further,  the  error  is  proportional  to  t/(£+k).  This  is  of 

+  +  + 

practical  interest  forthe  concentrations  of  0  ,  H  ,  and  N 
because  i  is  proportional  to  the  molecular  densities,  whereas 
k  is  proportional  to  atomic  densities.  Thus  the  error  in  Z^ 
is  of  order  (£/k)  compared  to  the  error  in  Y^,  i.e.  in  the  ratio 
of  molecular  to  atomic  densities.  Since  this  ratio  decreases 
rapidly  with  altitude,  the  ratio  of  the  error  in  Z ^  compared 
to  the  error  in  decreases  rapidly  with  altitude.  Finally, 
it  should  be  noted  that  transport  effects  are  needed  to 


reconcile  the  difference  between  the  observed  value  (X^) 


and  Z  ^ ,  not  the  difference  between  X^  and  Y^. 

It  is  useful  to  discuss  the  definition  of  chemical  life¬ 
times.  If  the  equation  of  continuity  does  not  involve  coupling 
to  other  ions,  the  chemical  lifetime  is  uniquely  defined  as 
the  reciprocal  of  the  factor  multiplying  the  ion  concentration 
in  the  loss  term.  When  ions  are  coupled  in  the  equation  of 
continuity,  there  are  two  lifetimes  that  can  be  defined.  In 
the  conventional  definition,  l/(£^+k^)  in  the  equation  for 
X |  is  the  lifetime  for  an  individual  ion.  To  discuss  the 
second  definition,  we  write  the  equation  for  X^  in  the  form 

°l  f  °2/U2+k2  +  t2)  *  "  klk2  /  («2+k2+t2)]  Xx  =  V-Fj  . 

It  can  be  argued  that  the  reciprocal  of  the  factor  multiplying 

X ^  in  this  expression  is  the  physically  meaningful  lifetime. 

In  the  absence  of  transport,  and  if  k2  <t2  r  the  lifetime 

then  becomes  l/^  ,  which  is  mucli  larger  than  l/(H^+kj)  in 

+  +  f 

the  upper  atmosphere  for  the  species  of  interest  (H  ,0  ,  H  ) . 
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3.  GENERAL  THEORY  OF  IONOSPHERIC  MOTIONS 


A.  Determination  of  Particle  Velocities 

The  drift  velocity  for  each  species  in  a  gas  mixture  is 
determined  by  the  momentum  transfer  equations  which  may  be 
written  as  (cf.  Schunk  1977) 


VPS  "  nsmsg  ' 


l-v 

n  e  (E+— v 
s  s  c  s 


x  B)  =  -Z  n  m  v  .  (v  -v,  ) 
t  s  s  st  s  t 


where  E  and  B  are  the  electric  and  magnetic  fields,  respectively, 
Ps,  ng,  mg,  eg  are  the  partial  pressure,  concentration,  mass, 
and  charge,  respectively,  g  is  the  gravitational  acceleration, 
and  v  is  the  collision  frequency  of  species  s  with  species  t. 
For  the  charged  species  it  is  convenient  to  deal  with  drift 
velocities  relative  to  the  neutral  medium: 


->■  t 

where  V  is  a  suitably  defined  velocity  for  the  neutrals.  Let  . 

F.  =  m.  g  +  e .  (E  +  — V  x  B)  -  Vp./n.  +  m.  1'  v.  (v  -  V) 
i  i  l  c  l  i  in  n 

n 

where  the  sum  over  n  indicates  that  the  summation  extends  over 
neutral  species  only.  Then 


Z  m.\>.  .(u.-u.)  +  Z  m.y.  u .  -  e.u.  x  B/c  =  ?. 
j  i  ID  i  3  n  1  in  1  11  1 

where  the  sum  over  j  extends  over  charged  species  only. 
Finally,  let 


v . 

i 


Z  v .  .  +  Z 
D 


n 


v . 

in 
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0>i  =  e^B/nu  C 
b  =  §/B 


where  B  is  the  magnitude  of  the  magnetic  field,  and  is  the 
gyrof requency.  Then 

miuA  -  mi  .5.  viA  +  miwi(b  x  «i>  =  ?i 

3^1  j  j 

To  solve  this  system  of  coupled  vector  equations,  define 
the  matrix  D  as  follows : 


.  v . 
1  1 


m 


~miV i j  (j^i) 


Then 

5  °i j  u'j +  mi“i(b  *  v  =  ?i 

The  solution  for  the  parallel  (to  B)  components)  follows 
immediately : 


b  •  u .  =  E  D..  b  •  F  . 

1  .  13  3 

1  J 

/v 

To  obtain  a  solution  for  the  transverse  components,  let  i2  be 

A  A  A  A 

an  arbitrary  unit  vector  normal  to  b,  and  i3  =  b  x  i2< 

A  /\ 

Equating  the  coefficients  of  i2 ,  i^  in  the  equation  of  motions. 

>:  D.  ^  u.„  -  m.u.u.  _  =  F .  _ 
j  13  J?  1  1  i3  i2 


>.  d:*  uj3  +  miMlu12  =  Pi3 
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These  equations  can  be  combined  to  give 

E  dJv  ~ —  orb  u.0  +  m.u>  .  u.  _  =  F.  „  +  E  D-.1  ~ —  F,  „ 

j  k  k  "V'k  k]  32  1  1  12  13  k  lk  Vk  k2 

Z(Z  d"*  ~ —  v7h  u  .  -  +  m.w.u.  .  =  -F.0  +  E  d'J  — —  F.  0 

j  k  k  Vk  kD  3  3  3  1  13  12  k  lk  Vk  k3 


Let  hT^  =  Z  D-3  —  d73  +  6.  .  m.w. 

3  k  lk  Vk  k3  13  i  i 


Then 


u.0  =  E  H .  .F. ,  +  E  P..  F. 

12  j  ID  i3  j  j2 

13  j  ID  32  j  id  33 


where  -  (S  Hifc  Dkj  > /nrw j 


The  total  transverse  component  is 


i0u._  +  i,u. ,  =  E  H . . (i_F . _  -  i,F._)  +  E  P..(i„F  +  i ,F  .  ,) 

2  i2  3  x 3  ^  i]'  2  j3  3  j2  ij  2  j2  3  j3 


b  x  (u-  x  b)  =  E  H. . (F .  x  b)  +  E  P. .  b  x  (F.  x  b) 
1  j  13  3  j  13  3 


=  (E  H.  .?.)  x  b  +  b  x  [(E  P...F.)  x  b] 


13  3 


13  3 


which  is  a  form  that  makes  no  eference  to  any  particular 

coordinate  system. 

,  ► 

The  definition  of  F^  , 


?i  s  V  +  eAc  -  Vp./ni  +  m.  E  v.n(^  -  V) 
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indicates  that  contributions  to  the  drift  velocity  come  from 

the  gravitational  acceleration  g,  the  total  electric  field 

->•  1 
e  =  E+  —  VxB  , 

c  ' 

gradients  in  the  partial  pressures,  and  contributions  arising 
from  differences  in  the  velocities  of  the  neutral  species.  The 
last-named  contributions  are  small  and  will  be  neglected 
hereinafter.  It  is  useful  to  distinguish  the  different 
contributions  to  the  velocity  vr  by  writing 

Ui  =  u^g)  +  u^e)  +  U^V) 

For  the  parallel  components: 

/S  ->•  ^ 

b'u^g)  =  (b*g)  E  D^nu 

b-u.  (c)  =  (b«  e)  ED.  .e  . 

1  j  J 

b« (v)  =  -E  Di . (Vp ./n  . )  •  b 
j  J  I  J 

For  the  transverse  components: 

A  -►  ~  A 
b  x  [u^(g)  x  b]  =  (g  x  b)  EH.  .m.  + 

j  13  3 

b  x  [u^(e)  xb]  =  (e  xb)  E  H..e.  + 

j  13  3 

b  x  [u.  (V)  x  b]  =  -E  H.  .  (Vp./n.)  x 
^  j  J  D  3 

B.  Current  Density 

The  current  density  J  is  given  by 


bx  (gxb)  E  P..m. 
*  ^  i]  D 


b  x  (e  x  b)  E  P^e. 


A  /S 


b  -  b  x  [E  Pi j (Vp j/n j )  x  b] 
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It  is  useful  to  write 


5  =  3(V)  +  3 (g)  +  3(e)  +  J(V) 


where,  owing  to  the  condition  of  local  neutrality. 


3iV)  =  V  -  E  n.e.  =  0 


c  i  i  1 


The  other  contributions  may  again  be  split  into  parallel  and 
transverse  components.  For  the  parallel  components: 


b  •  3(g)  =  (b-g)  i  E  n.e.  E  D. .m . 

i  j  J  J 

b  •  3(e)  =  (b*  e)  i  E  n.e.  E  D  .e. 

c  i  x  x  j  XJ  J 

i  ■  }(*)  =  I  I  Dij(VPj/nj)  •  £ 

For  the  transverse  components 

b  x  [3(g)  x  b]  =  (g  x  b)  i  E  n.e.  E  H. .m 

i  j  J  -* 

/v  ^  /s  i 

+  b  x  (gxb)— E  n.e.  E  P.  .m. 

ci  1  1  j  13  3 

b  x  [3(e)  x  b]  =  (e  x  b)  ^  E  n.e.  E  H. .e. 

c  i  x  j  J 

^  1 

+  b  x  (ex  b )— E  niei  E  P^^.e. 

°  i  j  J 

b  x  [3(V)  x  b]  =  i  E  n.e.  E  H.  .  (Vp./n.)  x  b 

c  ^  j  1  J  J  J 


The  conventional  electrical  conductivities,  which  deter¬ 
mine  the  contribution  to  the  current  density  from  electric 
fields,  are  seen  to  be 


Parallel:  cr^  =  —  E  n.e.  E  D.  .e. 

D  c  .  11.  IT  i 

1  3 


Hall: 


ot,  =  —  E  n .  e  .  E  H .  .e  . 


H 


c  i  1  1  j  ^  3 


Petersen:  a  =  —  E  n.e.  EP..e. 

p  c  .  11.  1]  ] 

i  3 

In  terms  of  the  electrical  conductivities,  we  may  write 


J(e)  =  o„(e-b)b  +  ou(e  x  b)  +  a  b  x  (e  x  b)  . 

D  ri  p 


C.  Ion  Drag 

The  force  exerted  on  the  neutral  constituents  by  charged 
particles  is  commonly  called  ion  drag.  Neglecting  differences 
in  the  neutral  particle  velocities,  we  may  write  the  force  (per 
unit  volume)  exerted  on  neutral  species  n  by  the  charged 
species  i  as 


?  . 

ni 


n  n  ni 


+  m  n  v  .  u. 
n  n  ni  i 


m. n . v  .  u. 
i  l  in  i 


by  the  symmetry  property  of  collision  frequencies. 

The  force  exerted  by  all  charged  particles  on  neutral 
species  n  is 
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?  =  Z  1  .  =  £  m. n.  v .  u. 

n  ini  ^  i  1  in  i 


Finally,  the  total  force  exerted  on  all  neutral  species  by 
all  charged  species  is 


?  =  £  In  =  £  n.m. (£  v.  )  u. 
„  •  i  i'  in'  i 

n  i  n 


Let  1  =  1(g)  +  1(e)  +  1(V) 


For  the  parallel  components: 


b*l  (g) 
b-?(e) 
b*l(V) 


(b'g)  £  n.m.  (Z  v.  )  £  D.  .m. 
*  .  ii  in  .  ill 


n 


(b.£)  z  n^llv  )  Z  D  e 

n  n  3 

A 

£  n.m.  (£  v.  )  £  D .  . (Vp  ./n . ) -b 
i  l  ivn  in  j  il  3  3 


For  the  transverse  components: 


b  x  [1(g)  x  b]  =  (b  x  g)  £  nim.  (£  vi  )  £  IT.m. 

i  n  j  J  J 


+  b  x  (g  x  b)  £  n.^nufZ  v^)  £  P.j.m 
i  n  j  J  3 


b  x  [1(e)  x  b]  = 


(b  x  e)  £  n.m.  (£  v  )  £  H.  .e. 

i  1  1  n  j  3  J 


+  b  x  (e  x  b)  £  n.^  (£  vi  )  £  P.  .e. 

i  n  j  3  J 

b  x  [?(V)  x  b]  =  £  n.m  (£  v  )  £  IT  (Vp./n.)  x  b 

i  n  j  -1  J  J 


+  £  n.m.  (£  v  )  £  P±.  b  x  [(Vp./n.)  x 

i  n  j  J  J 
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D.  Electric  Fields 


The  expression  relating  the  drift  velocity  of  a  charged 
species  to  the  electric  field  may  be  inverted  to  yield  an 
expression  for  the  electric  field.  For  the  parallel  component 

b-e  =  b*  u .  ( e  )  /  (E  D^-e.) 

j  J  J 

For  the  transverse  components: 

[  (E  H.  .e.)2  +  (E  P .  .e  . ) 2  ]  b  x  (e  xb)  = 

i  "i  i 

j  j  j  j 

=  -(E  H.  .e.)  [u.  (e)  x  b)  +  (E  P.  .e.)  b  x  [u.(e)  x  b] 

j1!!1  j1!!  -*■ 

A  vector  measurement  of  the  drift  velocity  ,  together 
with  measurements  required  for  the  calculation  of  ui(g)  and 
u^(V),  is  required  for  the  calculation  of 

Uf ( c )  =  ui  -  ui (g)  -  u^V)  , 

from  which  the  electric  field  can  be  computed.  This  procedure 
appears  to  be  feasible  in  regions  of  the  atmosphere  where 
one  ionic  species  dominates. 
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4.  THEORY:  ONE  IONIC  SPECIES 

Even  though  it  is  applicable  only  to  a  portion  of  the 
ionosphere,  the  theory  for  the  case  when  only  one  ionic  species 
and  electrons  are  present  helps  illustrate  several  features  of 
the  general  theory.  The  equations  of  motion  are 


m.  (v.  +  v .  ) 

x  in  le 

ui 

-  m. v . 
i  le  e 

A 

+  m^uK  (b 

x  u^ 

II 

►TJ4. 

H- 

m  ( v  +  v  .  ) 

e  en  ei 

-> 

ue 

-  m  v  •  u. 
e  ex  i 

A 

+  m  a)  (b 
e  e 

X  ue) 

II 

(0 

where  the  forces  F\  ,  ?  are  defined  by 

$  =  -Vp  /n+mg+ee,  s=i,e 

s  *s  s  s  s 

and  where  the  subscripts  i,  e,  and  n  denote  ions,  electrons, 
and  neutrals  respectively.  We  note  first  that  if  collisions  with 
neutrals  are  neglected  (v  =  =  0)  ,  then  the  equations  of 

motion  can  only  determine  the  velocity  difference  (u^  -  ue)  and 
not  the  velocities  of  each  species.  This  remains  true  when 
several  ionic  species  are  present. 


Parallel  Components 

It  is  useful  to  examine  the  solution  for  the  "uncoupled" 
case,  i.e.  when  ion-electron  collisions  are  neglected  (v^e  = 


b*  u.  -*•  (b-f* .  )/(m.v  .  )  h  V. 

i  l ' '  l  in  i 

b-ue  ->  (b-?e)/(meven)  -  Ve 
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Inasmuch  as  the  forces  and  Fe  are  of  the  same  order  of 
magnitude,  it  follows  that 

Vve  '  “  5  <Vre„,/(1Vin)  • 

The  "uncoupled"  case  also  serves  to  define  the  "uncoupled" 
diffusion  coefficient 


D„  =  kT  /m  v 
S  s'  s  sn 


In  the  presence  of  collisions,  the  solutions  are 


b-^  =  [(meven  +  mevei)  (b‘Fi)  +  (miVie)  (b‘Fe)  ]/D 

b-ue  =  Kmevei)  (b*?:i)  +  (mivin  +  mivie)  (^‘?e)  ]/D 


where 


D  -  ^mivin*  ‘Ven1  +  <Vin>  'Vei1  +  (mev.n>  <Vie> 


Since  n^  =  ng  in  the  present  case,  it  follows  from  the  symmetry 
property  of  collision  frequencies  that 


m.  v.  =  m  v  . 
1  le  e  ei 


Accordingly 


D  =  (m.  v.  )  (m  v  )  +  (m  v  .  )  [  (m. v .  )  +  (m  v  )  ] 
v  x  in'  '  e  en  e  ei  1  in'  e  en" 


=  (m.v.  )  [(mv  )  +  (mV  .  )  (1  +  a)  ) 

l  in  1  e  en  e  ei 
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and 


(b -  u±)  =  (mevei)  [b- +  Fe )  ] /D  +  (me^en)  (b-F±)/D 

(b-ue)  =  (mevei) [b* (#.  +  ?e)]/D  +  (mivin) (b-Fe)/D 

These  expressions  show  that  collisions  between  ions  and 
electrons  give  rise  to  a  velocity  component  that  is  common  to 
both.  It  should  be  noted  also  that  the  electric  field  does 
not  enter  this  component  because  it  appears  with  opposite 
signs  in  and  ?e.  The  "uncoupled"  velocities  ,  V are 
still  present,  but  they  are  divided  by  the  factor 


D/[(mivin)  (me%)]  =  1  +  (vei/ven)(1+a|* 
It  is  useful  to  note  that 


=  <Ven,/l,Vin> 


10 


-2 


at  all  altitudes.  The  parameter  a  is  a  measure  of  the  strength 
of  electron-neutral  coupling  compared  to  ion-neutral  coupling. 

A  second  useful  parameter  is 


v  •  /  [v  + 
ei  en 


( 1+a) v  •  ]  . 
ei 


As  the  fractional  ionization  (n^/n^  increases,  the  ratio 

(v  . /v _ )  increases  and  so  does  x.  It  is  convenient  to  refer 

ei  en 

to  x  as  the  (collis ional)  coupling  parameter.  The  reason  for 
this  choice  becomes  apparent  when  the  veloci  .ies  are  expressed 
as 

(b*ui)  =  (1-ax)  Vi  +  axVe  =  Vi  +  ax(Vg  -  V±) 

(b-ue)  =  (1-x)  Ve  +  xVi  =  Ve  -  x(Ve  -  Vi) 


l 
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where  and  Vg  are  the  uncoupled  velocities  previously 

defined.  The  collisional  coupling  changes  both  velocities,  the 

change  being  proportional  to  the  coupling  parameter  x  and  to 

the  difference  (V  -V^)  of  the  uncoupled  velocities.  The  change 

in  the  ion  velocity,  however,  is  smaller  than  the  change  in  the 

-2 

electron  velocity  by  a  factor  of  a(<10  ).  This,  of  course,  is 

the  result  of  the  fact  that  the  coupling  of  the  electrons  to 
neutrals  is  weaker  (by  a  factor  of  a)  than  the  coupling  of  ions 
to  neutrals.  In  fact,  as  the  coupling  parameter  approaches 
unity,  the  electron  velocity  tends  to  a  value  unrelated  to  its 
"uncoupled"  value. 

Since  collisions  between  particles  do  not  constitute 
external  forces,  the  total  force  per  unit  volume  remains  un¬ 
changed  by  collisions.  In  the  present  case,  this  force  is 


n. (b-F. )  +  n  (b • F  )  =  n.m. v.  (b*u. )  +  n  m v  (b*u  ) 
li  e  e  lx  in  1  e  e  en  e 


=  nimivin  [  (b-u^)  +  a (b • ug) ] 


=  n.m.v.  (V.  +  aV  ) 
1  1  in  i  e 


=n.m.v.  V.  +nmv  V 
l  i  in  i  e  e  en  e 


Thus,  irrespective  of  the  magnitude  and  nature  of  the 
collision  frequencies, 


(b*u.  )  4  a(b*u  )  =  V.  +  aV 
l  e  i  e 


In  some  theoretical  formulations  (see  for  example 

Schunk  and  Walker  1970)  it  is  assumed  that  (m  v  u  )  is 

e  en  e 
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of  order  (m  /m.  )  compared  to  (m.v.  u.)  ,  and  therefore  negli- 
6  j-  i  x  n  x 

gible.  This  assumption  is  not  correct  in  general.  The  ratio 
in  question  is  in  fact 


|<*V  -ax(V  -vi)  |/|  V±  +  ax(V  -V.)  | - >  f  otV  |  / 1  V±  | 

x+0 


x-*-l 


aVi|/|Vi+aVe| 


Since  | aV  | / | |  't  1  ,  as  shown  previously,  the  assumption  is 
incorrect  unless  the  collisional  coupling  parameter  x  is  near 
unity.  Data  from  daytime  observations  (listed  in  this  report), 
indicates  that  x  >  0.95  above  270  kilometers  and  x  >  0.99 
above  330  kilometers.  At  night,  the  fractional  ionization  is 
lower,  and  these  values  of  the  coupling  parameter  are  reached 
above  altitude  limits  which  are  'vlOO  kilometers  higher  than  in 
the  daytime. 

It  is  interesting  to  consider  the  effects  of  ion-electron 
collisions  on  the  parallel  component  of  the  current.  This 
is  proportional  to 


b- 


(Vj-Vfi)  ( 1-x-ax) 


-  <WX 


{v  /  [v 
en  en 


vei(1 


+  a)  ]  } 


As  x  approaches  unity  (which  is  equivalent  to  saying  as  electron- 

ion  collisions  dominate  electron-neutral  collisions)  the  effect 

of  electron-ion  collisions  is  to  reduce  the  parallel  conductivity 

by  a  factor  of  (ven/vei)  compared  to  the  "uncoupled"  case.  This 

result  further  shows  that  as  (v  /v  . )  tends  to  zero,  ions 

en  ei 

and  electrons  tend  to  move  with  the  same  velocity. 


i 
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Transverse  Components 


By  combining  the  equations  for  and  we  obtain 
uncoupled  equations  of  the  form 

^  A  ^  A  A  ^ 

A  u  +  A. (bxu  )  +  A_  b  x  (bxu  )  =  C  ,  s=i,e 
o  s  X  S  4  s  s 

where 


Aq  =  ur^(l  +  f  +  af )  ,  A^  =  -r(l-a),  A 2  =  -1 

C.  =  [  u  r  ( 1  +  f )  F .  +  arf  £  +  £ .  x  b  ]  /  ( m .  aj .  ) 

1  l  e  l  ’'li' 

Ce  =  [ r ( 1+af )  +  arf  E\  -  x  bJ/lnruK) 


and 


a 


=  (m  v  /m.v.  ),  r  =  (v.  /w.) 
e  on'  i  in  ’  m  i 


f  = 


(v  ./v  ) 

'  ei  en' 


The  parameter  a  has  been  discussed  before.  It  ranges  from 
3. 3x10  3  at  650  km  to  1.2xl0~3  at  160  km  in  the  data  used  in 
this  report.  The  ratio  r  of  the  ion-neutral  collision  frequency 
to  the  ion  gyrof requency  has  the  value  0.06  at  160  km  and 
decreases  with  altitude  essentially  in  proportion  to  the 
atmospheric  density.  Like  a,  r  is  insensitive  to  atmospheric 
conditions  where  its  value  is  significant.  The  parameter  f 
is  a  measure  of  the  fractional  ionization  and  ranges  from  0.5 

3 

at  160  km  to  4  x  10  at  657  km  in  our  data.  It  is,  of  course, 
highly  variable,  and  substantially  lower  at  night  than  in  the 
daytime.  Inasmuch  as  f  appears  in  a  product  with  r,  and  rf  is 
always  a  small  number  (>0.03),  knowledge  of  the  actual  value 
of  f  is  not  needed.  Taking  into  account  the  numerical 
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characteristics  of  the  parameters,  we  may  write  the  solutions 
in  a  simple  form  that  is  highly  accurate  at  all  altitudes: 


(1+r2)  (mia)i) 


bx(u^xb) 


[F.-ar2fF  ]  x  b  +  r  b  x  (?.xb) 

+  arf  b  x[(?.+?  )  x  b] 

1  e 


(1+r2)  (nua^)  b  x  (uexb)  =  -(1+r2)  (jr^xb)  -  ar2f[F’i+2Fe]  x  b 

+  ar  b  x  (Fexb)  +  arf  b  x  [(F^+F^Jxb] 


It  should  be  noted  that  the  "Hall"  components  (i.e.  the 

/s 

components  which  involve  a  single  vector  product  with  b)  are 

are  essentially  diagonal,  the  off-diagonal  elements  being  smaller 
2  -4 

by  a  factor  ar  f(<10  ).  The  diagonal  elements  are  essentially 

equal  to  l/(nuuK)  =  c/e B.  The  Hall  components  remain  diagonal 
even  when  several  ionic  species  are  present. 

On  the  other  hand  the  "Pedersen"  off-diagonal  components 
are  of  order  af (=vie/v^n)  compared  to  diagonal  components, 
and  this  product  may  be  small  or  large.  Finally,  the  magnitude 
of  the  Pedersen  components  is  of  order  r  (=v^n/u>^)  compared 
to  the  Hall  components. 

It  is  useful  to  examine  the  velocity  difference,  which  is 
proportional  to  the  current 

(1+r2)  (iikuk)  b  x  [(S.-Sj  x  b]  = 

(l+ar2f) (F. +?  )  x  b  +  r2  (?  xb)  +  r  b  x  [(?.-a?  )  x  b] 
l  e  e  l  e 

The  leading  "Hall"  component  is  independent  of  the  electric 

field,  since  and  ?  contain  the  electric  field  with  opposite 

signs.  Thus  the  Hall  electrical  conductivity  is  proportional 
2  2 

to  r  /(1+r  ),  and  decreases  very  rapidly  with  increasing 
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altitude.  On  the  other  hand,  the  Pedersen  electrical  con¬ 
ductivity  is  linearly  proportional  to  r,  and  decreases  less 
rapidly  with  altitude.  Finally,  we  note  for  convenience  that 


arf  = 


2  2,  ,  ,2. 
arf  =  (v.  v,  /to.  ) 
le  in  x 
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5.  RESULTS  AND  DISCUSSION 


In  this  section  we  present  some  results  using  the  formal¬ 
ism  of  Section  3  and  the  data  given  in  Appendix  A. 

A.  General 

It  was  shown  in  Section  3  that,  when  collisions  between 
charged  particles  are  taken  into  account,  the  velocity  of  a 
given  charged  particle  species  depends  on  the  forces  acting  on 
all  the  charged  species.  Figures  4,  5,  and  6  demonstrate  this 
for  the  parallel  component  of  the  velocity  of  0+  ions,  elec¬ 
trons,  and  N+  ions,  respectively.  More  specifically,  the  various 
curves  in  each  figure  represent  the  velocity  generated  in  the 
species  of  interest  by  a  force  on  another  species.  Before 
discussing  each  figure  separately,  we  note  that  these  figures 
give  the  velocity  (in  cm/ sec)  generated  by  a  force  of 
1.66  x  10  ^  dyne. 

Fig.  4  pertains  to  0+ ,  which  is  the  major  ion  through  most 
of  the  altitude  region  considered.  The  curve  labeled  "un¬ 
coupled"  represents  the  response  of  0+  ions  to  the  force  exerted 
on  them  in  the  absence  of  collisions  with  other  charged  parti¬ 
cles.  It  is  defined  explicitly  in  Section  A.  The  curve 
labeled  0+  represents  the  response  of  0+  ions  to  the  force 
exerted  on  them  in  the  presence  of  collisions  with  other 
charged  particles.  We  note,  as  in  Section  4,  that  the  intro¬ 
duction  of  collisions  has  a  very  small  effect  on  this  response 


function,  owing  to  the  much  higher  mobility  of  the  electrons. 

The  curve  labeled  "electrons"  represents  the  response  cf  0+ 

ions  to  the  force  exerted  on  electrons.  Because  0+  is  the 

major  ion,  this  response  function  is  nearly  identical  to  the 

response  function  to  the  force  exerted  on  0+.  Thus,  a  far 

as  0+  is  concerned,  the  presence  of  collisions  has  the  effect 

of  subjecting  0+  to  the  sum  of  the  forces  on  0+  and  electrons. 

The  effect  of  the  forces  on  the  other  (minor)  ions  is  an  order 

of  magnitude  smaller  in  the  altitude  range  considered.  This 

+  + 

result  remains  valid  even  below  220  km,  where  C>2  and  NO  have 
higher  concentrations  than  0+,  because  at  these  altitudes 
collisions  with  neutral  particles  play  a  major  part,  and  the 
ions  are  not  strongly  coupled. 

Fig.  5  shows  the  response  of  electrons  to  forces  exerted 
on  each  charged  species.  As  pointed  out  in  Section  4,  the 
introduction  of  collisions  reduces  drastically  the  response 
of  electrons  to  the  force  exerted  on  them.  This  is  evident 
in  Fig.  5  from  a  comparison  of  the  curves  labeled  "electrons" 
and  "uncoupled"  (no  collisions) .  With  increasing  altitude, 
collisions  with  0+  predominate  over  collisions  with  neutrals, 
and  the  motion  of  electrons  is  effectively  governed  by  the  sum 
of  the  forces  on  0+  and  on  electrons.  The  effect  of  minor  ions 
is  an  order  of  magnitude  smaller. 

Fig.  6  represents  the  respone  of  N+  ions  to  forces  exerted 
on  each  charged  species.  N+  is  a  minor  ion  throughout  the 
altitude  range  considered.  Fig.  6  shows  that  the  motion  of  N+ 
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ions  is  controlled  by  the  sum  of  the  forces  on  0+  (the  major 
ion)  and  electrons.  The  response  to  the  force  on  N+  is 
negligible  compared  to  its  value  in  the  absence  of  collisions 
(as  depicted  by  the  curve  labeled  "uncoupled"),  except  at  the 
lowest  altitudes,  where  collisions  with  neutrals  dominate. 

At  the  upper  altitudes,  where  collisions  with  0+  are  dominant, 
the  response  functions  of  N+  are  identical  to  the  response 
functions  of  0+,  which  means  that  the  motion  of  N+  ions  is 
identical  to  the  motion  of  0+  ions.  This  result  is  valid  for 
all  the  minor  ions. 

We  can  summarize  the  results  for  the  parallel  components 
as  follows: 

-  Where  charged  particle  collisions  dominate  ion-neutral 
collisions,  minor  ions  move  with  the  same  velocity  as  the  major 
ion; 

-  Except  as  noted  below,  the  velocity  of  the  major  ion  and 
of  electrons  is  proportional  to  the  sum  of  the  forces  exerted 
on  the  major  ion  and  on  electrons; 

-  In  the  presence  of  parallel  electric  fields,  the  velocity 
of  electrons  may  be  substantially  different  from  the  velocity 

of  the  ions.  The  effect  of  electric  fields  on  ions  is  of 

-2 

order  a(=m  v  /m.v.  <10  )  compared  to  the  effect  of  electric 

e  en  1  in 

fields  on  electrons. 

For  the  transverse  velocity  components  we  note  that  the 
"Hall"  contribution  is  negligible  except  for  a  force  due  to  a 
transverse  electric  field,  and  that  the  "Pedersen"  contribution 
is  negligible  in  the  altitude  range  considered. 


B.  Calculation  of  0+  Drift  Velocity. 

We  have  calculated  the  drift  velocity  of  0+  ions  using  the 
equations  of  Section  3.  In  the  expression  for  the  force  on 
species  s 

?  =mg+eE-Vp/n 
s  s3  s  rs  s 


we  have  assumed  that  the  gradient  of  the  partial  pressure  has 
a  vertical  component  only.  Further,  in  the  absence  of  any 
data  on  electric  fields,  we  have  calculated  separately  the 
contribution  from  parallel  (to  the  magnetic  field)  and  eastward 
electric  fields.  The  results  are  shown  in  Figures  7,  8,  and  9. 

Figure  7  shows  the  effect  of  the  gravitational  force  on 
charged  particles,  and  demonstrates  that,  in  the  region  where 
collisions  between  charged  particles  dominate  collisions  with 
neutrals,  the  gravitational  force  affects  all  species  equally. 

Figure  8  shows  the  contribution  to  the  velocity  of  0+  from 
various  terms.  The  vertical  component  of  the  velocity  is  given 
by  sin  I  times  the  value  shown  in  the  figure  (where  I  is  the 
magnetic  inclination) .  Inasmuch  as  the  magnitude  of  the  velocity 
varies  by  104  over  the  altitude  range,  a  logarithmic  plot  has 
been  used,  and  a  dashed  curve  indicates  negative  contributions. 
The  curve  labeled  -G  represents  the  contribution  from  the  gra¬ 
vitational  force,  as  in  Fig.  7.  The  curve  labeled  N  (or  -N) 
represents  the  contributions  from  gradients  of  the  concentra¬ 
tions  of  charged  particles.  The  curve  labeled  T  (or  -T) 
represents  contributions  from  the  gradients  of  electron  and 


ion  temperatures.  All  ions  were  assumed  to  be  at  the  same 
temperature.  The  curve  labeled  S  represents  the  sum  of  the 
aforementioned  contributions.  Figure  8  indicates  that  the 
contribution  from  temperature  gradients  is  relatively  small, 
that  the  contribution  from  the  gravitational  force  is  the 
largest  contribution  over  most  of  the  altitude  ranges  and 
that  the  sum  of  the  contributions  is  negative  below  about  590  km. 

Fig.  9  displays  various  contributions  to  the  vertical  com¬ 
ponent  of  the  0+  drift  velocity.  As  in  Fig.  8,  a  logarithmic 
plot  has  been  used,  and  dashedlines  indicate  negative  values. 

The  curve  marked  ug  represents  the  contributions  to  the 

vertical  velocity  from  the  gravitational  force  and  from  the 

2 

gradients  of  partial  pressures.  u_  is  proportional  to  sin  I, 
where  I  is  the  magnetic  inclination,  and  thus  tends  to  zero 
at  about  375  km,  which  occurs  above  the  magnetic  equator. 

Also  shown  in  Fig.  9  by  the  curve  labeled  is  the  contri¬ 

bution  to  the  vertical  velocity  of  a  (constant)  parallel  elec¬ 
tric  field  of  1  millivolt  per  meter.  This  contribution  is 
proportional  to  sinl,  and  also  vanishes  at  the  magnetic 
equator.  The  curve  labeled  represents  the  "Hall"  contri¬ 
bution  to  the  vertical  velocity  of  a  geomagnetical ly  eastward 
electric  field  of  3  millivolts  per  meter.  It  should  be  noted 
that,  insofar  as  the  calculation  of  drift  velocities  is  con¬ 
cerned,  this  contribution  is  the  same  for  all  charged  particles, 
including  electrons.  Finally,  the  filled  circles  represent  the 
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vertical  component  of  plasma  drift,  as  measured  by  the  RPA 

experiment  on  the  AE-C  satellite  (Hanson,  1973) .  A  comparison 

of  the  measurements  with  the  various  curves  indicates  that  an 

eastward  electric  field  of  between  3  and  5  millivolts  per  meter 

is  sufficient  by  itself  to  explain  the  observations.  Further, 

it  appears  that  if  a  parallel  electric  field  is  present,  it 

must  be  less  than  ^0.01  millivolts  per  meter  over  most  of  the 

altitude  range.  Up  to  about  500  km,  the  contribution  from 

dominates  ug,  and  the  calculated  drift  velocity  can  be  said 

to  be  in  harmony  with  observations.  Above  that  altitude, 

however,  the  values  of  u  dominate.  It  should  be  remembered 

s 

that,  in  calculating  gradients,  it  was  assumed  of  necessity 
that  variations  occur  only  in  the  vertical  direction,  and 
that  no  ( geomagnetical ly)  north-south  variations  are  present. 

It  may  very  well  be  that  this  is  not  a  valid  assumption  in 
the  equatorial  region. 

C.  The  Equation  of  Continuity  for  0+ 

To  provide  a  satisfactory  explanation  of  the  observed 
altitude  profile  of  0+  ions,  we  must  be  able  to  calculate  the 
divergence  of  the  flux  in  the  equation  of  continuity 

P  -  L  =  V  •  (nu) , 

in  regions  where  transport  is  significant.  This  task  requires 
knowledge  of  the  variation  of  the  concentration,  and  of  the 
drift  velocity,  in  a  three-dimensional  frame.  Unfortunately, 
the  data  do  not  allow  the  calculation  of  horizontal  derivatives. 
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and  the  assumption  that  only  vertical  variations  are  present 
is  not  consistent  with  the  equation  of  continuity. 

D.  Conclusions 

The  motion  of  charged  particles  in  the  upper  atmosphere 
has  been  investigated  with  the  aid  of  satellite  data.  It  has 
been  shown  that  the  major  ion,  0+,  determines  the  motion  of 
all  minor  ions.  The  observed  plasma  drift  velocity  in  the 
vertical  direction  is  in  harmony  with  the  assumption  of  an 
eastward  electric  field  of  between  3  and  5  millivolts  per 
meter.  It  appears  that  the  assumption  that  horizontal  varia¬ 
tions  are  negligible  is  not  valid  in  the  present  instance, 
owing  probably  to  the  fact  that  the  altitude  region  where 
transport  is  important  lies  within  15°  of  the  magnetic  equator. 
The  present  investigation  clearly  demonstrates  the  need  for 
electric  field  measurements  and  measurements  of  plasma  drifts. 
The  approach  outlined  in  this  investigation  will  be  applied  to 
midlatitude  data,  where  the  assumption  that  horizontal  varia¬ 
tions  are  negligible  is  valid. 


6.  ELECTRON  VELOCITY  DISTRIBUTION 

Photoelectrons  and  secondary  auroral  electrons  play 
an  important  role  in  the  coupling  of  the  solar  ionizing 
ultraviolet  flux  and  the  fast  primary  particles  to  the 
upper  atmosphere.  Theoretical  studies  of  the  transport 
and  thermalization  of  the  photoelectrons  and  secondary 
auroral  electrons  are  necessary  to  determine  how  the  incident 
energy  gets  partitioned  among  heat,  ionization,  dissociation, 
luminosity  and  other  modes  in  the  atmosphere.  The  presence 
of  a  parallel  electric  field  may  influence  the  partitioning 
of  the  energy  and  therefore  provide  the  possibility  of 
developing  diagnostics.  Our  preliminary  effort  to  study 
the  effects  of  electric  fields  on  the  electron  velocity 
distribution  has  been  to  modify  and  generalize  an  accurate 
electron  deposition  code  (Victor,  Kirby-Docken  and  Dalgarno 
1976) .  with  these  modifications  we  had  hoped  to  develop 
a  perturbation-iterative  method  to  introduce  the  effects 
of  weak  electronic  fields.  Limits  of  time  and  funding 
have  limited  progress.  The  code  takes  explicit  account 
of  the  discrete  nature  of  the  electron  energy  loss  process 
and  employs  cross  section  data  based  on  a  recent  critical 
reivew.  Local  deposition  is  assumed  so  that  accurate 
results  can  be  obtained  only  for  altitudes  below  about  300  km 

The  latest  original  version  of  the  code  was  operational 
on  a  )’.9  computer  system  at  NASA  Goddard  Space  Flight  Center. 
Minor  modifications  were  made  so  that  the  code  would  run 
efficiently  on  CDC  6600  series  computers  such  as  those 
at  AFGL. 
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The  code  was  then  generalized  so  that  it  would  calculate 
the  steady  state  velocity  distribution  for  arbitrary  initial 
velocity  (E  _<  200  eV)  distributions  obtained  from  a  file 
(TAPE4)  read  by  subroutine  SETPRO.  A  card  copy  of  the  source 
deck  and  the  cross  section  data  was  supplied  to  AFGL.  A 
listing  of  this  deck  is  given  in  Appendix  C.  A  listing  of  a 
typical  production  file  (TAPE4)  is  given  in  Appendix  D.  A 
small  code  to  produce  an  input  file  appropriate  to  an  auroral 
secondary  electron  production  distribution  (using  the  energy 
distribution  given  by  Opal,  Peterson  and  Beaty  1971)  has  been 
written.  A  listing  of  the  source  code  and  sample  output  is 
shown  in  Appendix  E. 
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ALTITUDE  (km) 


l.  The  concentration  of  molecular  ions  plotted  aqainst 
altitude.  Dots  indicate  measurements  by  the  MIMS 
experiment  (Hoffman  1973)  on  orbit  594  (uplot)  of  the 
AE-c  satellite.  Pluses  indicate  theoretical  values. 
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text  for  units) .  The  curve  labeled  "uncoupled  reters  to 
to  the  force  exerted  on  them  in  the  absence  of  collisions. 
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log  (velocity  in  cm /sec) 

Fig.  7.  Contribution  to  the  drift  velocity  of  ions  and  electrons  from  ,the  gravitational 
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Fig.  9.  The  vertical  component  of  the  0+  drift  velocity.  u?  represents 
contributions  from  the  gravitational  force  and  gradients  of 
partial  pressures.  En  and  E*  represent  contributions  from  a 
parallel  and  eastward  electric  field  of  1  and  3  mV/meter, 
respectively.  The  filled  circles  represent  the  observed  vertical 
component  of  plasma  drifts. 
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APPENDIX  A 


DATA 

The  data  utilized  in  this  report  pertain  to  the  upleg  of  orbit 
594  of  the  Atmosphere  Explorer  C  satellite  (Dalgarno  et  al. , 

1973) .  The  measurements  made  on  this  occasion  (February  8, 

1974)  are  representative  of  an  undisturbed  daytime  atmosphere  in 
the  mid-latitude  and  equatorial  region.  The  data  are  summarized 
in  three  tables. 

Table  A-l  lists  ion  concentrations  as  a  function  of  altitude 

Z  (in  km)  as  measured  by  the  MIMS  instrument  (Johnson,  1973) . 

-  3 

Concentrations  are  in  units  of  cm  and  listed  in  logarithmic 
(base  10)  form.  The  electron  concentration  is  the  sum  of  the 
ion  concentrations.  An  entry  of  -5  signifies  that  no  measurement 
of  the  concentration  was  made  at  that  altitude. 

Table  A-2  lists  neutral  concentrations  (in  cm  )  in  loga¬ 
rithmic  form.  For  N2,  He,  0,  and  N(^S)  the  measurements  by  the 
OSS  instrument  (Nier,  1973)  were  used  to  establish  a  Bates-Walker 
(Walker,  1965)  profile  from  which  values  were  obtained  at  all 
altitudes.  The  photochemical  model  of  Oppenheimer  et  a]L.  (1977) 
was  used  in  conjunction  with  a  Bates-Walker  profile  to  establish 
the  O2  concentration.  The  H  concentration  represents  a  model 
calculation  (Hedin  et  a_l.  ,  1977)  .  The  NO  concentration  was 

measured  by  the  UVNO  experiment  (Barth  et  al_.  ,  1973)  and  was 

2 

extrapolated  to  altitudes  above  260  km.  The  N(  D)  concentration 
was  measured  by  the  VAE  experiment  (Hayes  et  al_ . ,  1973)  and 
was  extrapolated  above  300  km. 


Table  A-3  lists  the  electron  temperature  (TE)  measured  by  the 
CEP  experiment  (Hanson  et  al.  ,  1973);  the  ion  temperature  (TI) 
measured  by  the  RPA  experiment  (Brace  et  al. ,  1973),  the  neutral 
temperature  (TBW)  inferred  from  the  0  and  ^  concentration  pro¬ 
files,  the  zenith  angle  (SZA) ,  magnetic  latitude  (MLAT) ,  latitude 
(LAT) ,  and  local  solar  time  (LST) .  No  measurement  of  the  neutral 
temperature  (TN)  was  available  on  this  occation.  Temperatures 
are  given  in  degrees  Kelvin,  angles  in  degrees,  and  time  in 
ho  urs . 
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TABLE  A-l.  Logarithmic  Concentrations  of  Ions  and  Electrons  (cm-3) 
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TABLE  A-2 .  Logarithmic  Concentrations  of  Neutrals  (cm  ) 
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TABLE  A-3.  Temperatures  and  Orbital  Parameters 
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APPENDIX  B 


COLLISION  FREQUENCIES 

This  appendix  contains  expressions  for  the  calculation  of 
various  types  of  momentum  transfer  collision  frequencies.  These 
are  defined  in  such  a  way  that 


-"s'Vst'V'V 


represents  the  force  per  unit  volume  exerted  by  species  t  on 
species  s.  n  ,  m  ,  and  u  are  the  concentration,  mass,  and 
velocity,  respectively,  for  species  s,  and  v  t  is  the  momentum 
transfer  collision  frequency  of  species  s  with  species  t.  As 
defined,  the  collision  frequencies  satisfy  the  symmetry  property 


n  m  v  =  n  m  v 

s  s  st  t  t  ts 


The  preceding  definition  is  the  most  frequently  used  (cf. 
Schunk,  1977).  It  should  be  noted,  however,  that  other  defini¬ 
tions  are  in  use  as  well.  For  example.  Banks  and  Kockarts 
(19  7  3)  ,  who  have  given  a  comprehensive  review  of  collision 
frequencies,  define  them  in  such  a  way  that 


v  (Banks)  =  (mg/ugt)  vgt(Schunk)  . 
where  p  =  m  m  /(m  +m  )  is  the  reduced  mass. 

St  St  St 


Electron-Neutral  Momentum  Transfer  Collision  Frequencies 
v(e-O)  =  8.2  x  10-10  n  (0) 

v.(e-N2)  =  2.33  x  10-11  n  (N2)  ( 1- 1 . 21xl0_4Te)  Tq 


v(e-Oz)  =  1.82  x  10  xu  n(02) (1  +  3.6  x  10  T  )  Te 
v(e-He)  =  4.6  x  10_1°  n(He)  Te1/2 

v(e-H)  =  4.5  x  10-9  n(H)  (1-1.35  x  10_4  Tg)  Te1/2 
v(e-N)  =  6.0  x  10_1°  n(N)  Te1/2 

where  T^  is  the  electron  temperature,  n(x)  is  the  concentration 
of  species  x  in  units  of  cm  ,  and  the  collision  frequencies 
are  in  sec  .  Except  for  the  collision  frequency  with  N,  which 
is  an  interpolation  based  on  the  values  for  He  and  O,  these 
expressions  are  given  by  Banks  and  Kockarts  (197  3)  and  they  per¬ 
tain  to  elastic  collisions.  These  authors  estimate  that  the 
accuracy  for  N2  and  0 2  is  about  20%,  for  H  25%,  for  He  10%,  and 
state  that  the  result  for  0  is  uncertain  both  as  to  magnitude 
and  as  to  the  dependence  on  Tg. 

Ion-Neutral  Momentum  Transfer  Collision  Frequencies 

At  low  temperatures  the  ion-neutral  interaction  arises  from 
the  induced  dipole  polarization  of  the  neutral  by  the  ion.  The 
corresponding  polarization  collision  frequency  is 

v.  =  2.6  x  10~9  n (n)  (A.  a  ) 1/2/A. 

in  m  n  i 

_  3 

where  n(n)  is  the  neutral  concentration  in  cm  ,  A^n  =  A^A^/fA^+A^ 

is  the  reduced  mass,  A.  and  A  are  the  ion  and  neutral  masses  in 

atomic  units,  a  is  the  neutral  polarizability  in  units  of 
-24  3  -1 

10  cm  and  v^n  is  in  sec  .  Except  for  the  ad3ustment  re¬ 
quired  by  the  different  definitions,  these  expressions,  as  well 
as  the  polarizabilities  listed  below,  are  given  by  Banks  and 


Kocka  rts  (  197  3)  . 


Polarizabilities  of  Neutral  Gases 


Neutral  Gas 
H 

He 

N 

0 

°2 

NO 


-24  3 

a  (10 _ cm  ) 

0.667 

0.82 

1.1 

0.79 

1.59 

1.74 

1.76 


At  higher  temperatures,  a  repulsive  force  counters  the 
polarization  force.  The  resulting  collision  frequencies  depend 
on  the  details  of  the  particular  ion-neutral  pair,  but  there  is 
no  data  for  collisions  at  high  temperatures.  This  lack  of  data 
is  compensated  to  some  extent  by  the  fact  that  in  the  high- 
temperature  regions  of  the  atmosphere  the  most  important 
collisions  are  resonant  charge  transfer  collisions,  i.e.  colli¬ 
sions  between  an  ion  and  its  parent  neutral.  The  values  given 
by  Banks  and  Kockarts  (1973),  adjusted  by  the  factor  p^/nn,  are 


v(0+-0)  =  2.4  x  10-13 

T1/2 

(10.6-0.67 

1O910 

T)  2 

n  (0 ) 

T>  470 

v (02-02 )  =  1-7  x  10"13 

T1/2 

(10.6-0.76 

1O910 

T)  2 

n(02) 

T>1600 

v(N2~H2)  =  1. 85  x 10-13 

Tv2 

(14.3-0.96 

1O’l0 

T)  2 

n(N2) 

T>  340 

v(H+-H)  =  9.5  x  10“13 

T1/2 

(14.4-1.17 

1O910 

T)  2 

n  (H) 

T>  100 

v (N+-N)  =  2.6  x  10~13 

t1/2 

(10.4-0.64 

1O910 

T)  2 

n  (N) 

T'  550 

v (He-N)  =  4. 85  x 10"13 

rpl/2 

(11.6-1.05 

log10 

T)  2 

n  (He) 

T''  100 
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where  T  =  is  the  sum  of  the  ion  and  neutral  temperatures. 

These  expressions  are  valid  for  temperatures  above  the  limit 
indicated  for  each  reaction. 


Coulomb  Collision  Frequencies 

For  charged  particles  the  collision  frequency  arising  from 
the  Coulomb  interaction  is 


“st  -  8-47  11  10'2  A  ’sW/VVsf2 


,  ~  2_  2,  1/2  3/2 

*  l-27  *s  zt  Ast  Wst 


with  £nA  =  15  (cf.  Banks  &  Kockarts,  1973).  nfc  is  the 
concentration  of  species  tf  A  ,  T  ,  and  z  are  the  mass  (in 

S  p  s 

atomic  units) ,  temperature  (in  °K)  and  charge  (in  atomic  units) 
of  species  s,  and 


Collision  frequencies  for  electrons  and  for  0+  are  shown  in 
Figure  B-l,  and  Figure  B-2,  respectively,  based  on  the  data 
given  in  Appendix  A. 
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Fig.  B-l.  Collison  frequencies  of  electrons  with  ions  and  neutrals 


(in  sec  .  Also  shown  is  the  sum  of  the  collision 
frequency  with  all  ions. 
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log  (collision  frequency) 
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frequencies  with  all  other  ions. 
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APPENDIX  C 

LISTING  OF  THE  ELECTRON  ENERGY  DEPOSITION 
CODE  AND  THE  CROSS  SECTION  DATA  FILE 


r>  n  n  o  n  n  n  r>  nr\n 


PROGRAM  RCTEMPR  (INPUT (OUTPUT ,TAPil ,TAPE2»CUTFUT,TAPE3»INPUT, 

1  TaPE4,TAP15, TAPES, TAPE7) 

****  RCTENPR  ****  770912 

THIS  IS  THF  DRIVER  THAT  IS  WRITTEN  TO  INTERFACE  WITH  THE 
PRIMARY  PRODUCTION  FIIES 
COMM0N/XSX5/  NOSPEC ,  NCEXC(4),  NGR1C(4,20), 

1EGRID(4,20,25) ,5PLC0(4, 20, 25) ,SPLC 1(4,20,25 ) ,SPLC2(4, 20 ,25), 
2SPLC3(4,20,25) , THRESH ( 4 , 20 } ,  MONFS  (  4  >  ,  MONGD  ( 4 ) , t 10NGD ( 4 , 25 ) , 
3DIFC0N(4),T1XS0(4,25),TIXSI (4 , 25 ) , T1 XS2 ( 4 , 25) , TIXS3< 4 , 25) ,C1XS0  l  4 , 
42b)  ,CIXS1 (4,25) ,D1XS2(4,25) ,CIXS3(4,25) ,  1HRJCZ(4#6) , 

5P  h ACT  (4,6,2b) 

COMMGN/UENS/TE,NOALT,IALT( tu) , XUENS(4) ,  xeden 
CQMMON/GRID/lNTGRin,XGRID(400)  , VGR ID ( 400 ) , FRGGR ID ( 400 ) 

COMMON  /PCHNT/IPCINT(4,20)  ,IIP0tiT(4)  ,  It  PONT  (4) 

C0*mCn/ALPHA/NAmE(4)  ,S1ATE{<I,20)  ,  1ST  AT EC  4, 6) 

CO  A  WON  IbGN,10(PL,IAVGF,FLXMULT 
HLw ARE ,  THE  OHDEfi  CF  SPECIES  IN  HECXSN  RUST  NATCH  WITH 
THOSE  IN  THE  PROCOCTION  FILE 
CALL  SFTCRID 

SETGRIO  SETS  UP  THE  INTEGRATION  GRID 
WHI IE(2,4) 

4  f ORMAT  ( *  ENTER  NUhPER  CE  ALTITUDE  VALUES! 12)*) 

REAU(3.1)  NOALT 

NOALT  IS  THfc  HUMBER  OF  ALTITUDE  VALUES  TC  BE  DONE 

1  FURR AT  112) 

*R  I TE ( 2 , 5 ) 

5  FORNAT  (  *  ENTER  LUMBERS  FOR  ALlITUOt  RECORDS  (  20  14  )*  ) 

READ (3. 2)  ( 1  ALT 1 1 ) , I* l, NOALl) 

wRIIE(2,b) 

6  FORWATC*  ENTER  n  U  m  H  E  R  CF  INTERVALS  FOR  AVG.U2)*) 

RE AD ( 3,1)  IAVGF 

N.B.  1AVGE.LT. 2  IwFUES  AVERAGED  FLUX  AND  FLUX*4P1  NCT  CCWFUTED 
WRITE (2, 10) 

10  EUPN AT ( *  ENTER  ELUX  W U LT I PL IC AT  ION  F ACTOR ( F 3 . 1 ) * ) 

REACl3.ll)  PLXMULT 

11  FORMAT! F 3.1) 

2  FORMAT (2014) 

CALL  REDXSN 
CALL  SETPRCC-l ) 

CO  3  I *1 , NT AL T 
CALL  SETPRC ( I  ALT ( I  )  ) 

CALL  WORK 

3  CONTINUE 
ENO 

SUBROUTINE  SETPRC (KEY  j 

COMMON  1BGN .TOTEI . IAVGP , ELXWULT 

C0NWCN/GRID/INTGPIC,XGRID(400)  ,  VGR1D1400)  .  FHCGR.TC  (  400  ) 

DIMENSION  XPRC1400) , YPRG(40u) ,C(400) 

DIMENSION  *aRCS(19) 

CORWOn/UENS/TE.NCALT, IALT( 10) ,X0tN5(4) .XEDEN 

THIS  ROLTI.'.E  SETS  UP  THE  PRIMARY  ELECTRON  PRODUCTION  DATA 
R  F Y  *•  1  TO  HEAD  PRODUCTION  GRID  DAT  A  ,  KEYNC  READS  THE  PPOCLCTION 
OATA  FOR  The  KtrtH  ALTITUDE  Cn  THE  FILE 
IF  (  KEY  ,'NE' ,  “  1  )  GO  TO  1 
READ  (4,92)  WORDS 
52  FORMAT  ( 1 9 A 4 ) 

WRITE ( 2 , 2)  wCRDS 


2  FURMAT  (1H1.19A4) 
READ  (4,52)  wORUS 
RE AO (4,  1)  NPTSPHC 
wHITF(2,2J)  NPTSPRC 
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23  F0RMAT(14,=  *  NP7SPR0*) 

3  FORMATt 15) 

PEAD(4,4)  (XPWO(I),  1*1,NP15PR0) 

WR1IE(2,21> 

21  FORMAT { *  ELECTRON  PRODUCTION  GRID  *) 

WRITE (2 ,22)  (XPRO(I),l*l .NPTSPRO) 

22  FORMAT  (IOE12.3) 

nmax«nptsprq*i 

DO  17  1*1 ,  NPT  5PRC 
17  XGRIO(I)sXFRO(l) 

CO  117  1*1,  iv MAX 

117  VGPID(I )=5.93E+7*SQHT(0.5*(XGRID(I)+XGRID(1+1) )  ) 

4  FORM AT ( 5E1 5  «  5 ) 

REWIND  4 
RETURN 

1  CONTINUE 
I*NPTSPR0/5 
F*  I 

F1=FLOAT(NPTSRRO)/5.0 
IE(F .EO.Fl )  GO  TC  5 
1*1  +  1 

5  CONTINUE 
1  =  1  +  3 

CO  6  J* 1 , I 
READ (4, 3) 
b  CONTINUE 

IF(KET.EO.l)  GO  TO  7 
IsNPTSPKO/t) 

F  =  I 

E  l=FLO AT  (NPTSPRO)/ 6,0 
lF(F.FC.Fl)  GO  TC  8 
1  =  1  +  1 

R  CONTINUE 

C  S1ATFMENT  1  =  1  +  10  REPLACED  8Y  1  =  1  +  11  (PC , 1 0/ l 7/75  ) 

1  =  1  +  1  1 
IL=KEY-1 
CO  9  0*1, 1L 
DO  9  E  =  1  , 1 
READU.3) 

9  CONTINUE 
7  CONTINUE 
READ(4, 10JZA 

10  FQRWAT(25X,E 12.b) 

RtA0(4,10)  ALT 

13  FORmaK*  ALTITUDE  iS*,Fl2.b) 

READ(4,10)  T 

WRITE  (5,13)  A  L  1 
WRITF  (2,11)  ALT»ZA,T 

11  FORMAT  (IX, 119(1 H“)/24H  ALT  SZA  TEMP/3F8.2) 

READ(4,14)  XCENSC4) 

RE  AC (4, 14)  X  C  t NS ( 2 ) 

READ (4,14)  XCENS(J) 

RE  AD ( 4 , 14  )  XIEN5(1) 

14  FORMAT(25X,Elb.b) 

READ (4, 15)  Tt 

15  FORMAT  (2bX,F16.b) 

READ (4, 20)  X  EDEN 

20  FORMAT  (22X,£lb.b) 

RE  AD ( 4 , 3 ) 

C  ANOTHER  DUMMY  READ  ECK  “TOTAL  ION  DENSITY*  ( RC , 1 0/ 1 7/75 ) 

READ  (4,3) 
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REAOU,  16)  ( YPRO(K)  ,K»1  , NPTSPRO) 

16  FORMAT  (6612,4) 

TOTEtsO.O 
DO  32  K«1  ,  N  M  A  X 
YPRO(K)wYPRC(K)*ELXMULT 
TOTEL»TOTEL*YPRO(K) 

PROGRIO(K)*YFRO(K)/(XGR1D(K*1 )-XGRlD(K)) 

C(K)»PROGRID(K)*VCRIO(K) 

C(K)»C(K)/12.568 
32  CONTINUE 

MR JTE ( 2 , 400  ) 

400  FORMAT!*  INITIAL  FLUX  *) 

WRITE (2,22)  ( C ( K ) , Ks 1 , NM AX  ) 

REWIND  4 

FIND  END  CF  PRODUCTION  SPECTRUM 
KxNPTSPRO-1 

18  IF(AB5(YPRClK)-0.0).GT.1.0E-6)  GO  TO  19 
K»K-1 

GO  TO  18 

19  CONTINUE 
IBGN=K 

WRI1E(2,31)  IBGN 
31  E0RMAT(I4,*  *  IBGN*) 

WRITE (7, 200)  ALT,T,XDFNS,TE.,XEDEN 

200  F0RMAT18E10.4) 

WHITE  (6.201)  ALT 
RETURN 

201  FORMAT  (F7.1) 

END 

SUBROUTINE  *ORK 

DIMENSION  DFPROS( 10) ,CFE(4.20) ,CFl(4,b) ,110(4) ,PRDFA(400) .PRCFASd 
10),FELS(10) ,EtL(400),PROBt(4,20),PPOBI(4,6),EEF(400),FELWP(400)  ,PR 
2C(400)  ,DilMMY(b)  ,0(13)  ,GE(20) 

COMMON  IBGN .TOTE (, IAVGE 

COMMGN/XSXS/  NOSPe;C,  nCEXC(4),  NGR1C  (4,20), 

1FGRI0(  4,20,25)  ,  SPl.CO  (  4 , 20 , 25  ) ,  SPLCl  (  4 , 20 , 25  )  ,  SPLC2  (  4 , 20 , 25  )  , 

2SPLC 3 (4, 20, 25) .THRESH (4, 20) , N IONFS ( 4 ) , N 10NGD ( 4 ) , EICNGD ( 4 , 25 ) , 

3D IFCON ( 4) ,T 1X50(4,25 ) ,11 XS1 ( 4 , 25 ) , T1XS2 ( 4 , 25) , TIXS3 ( 4 , 25 ) , C IXSO ( 4  , 
425),CIXSl(4,2b),DIXS2l4,25),DlXSJ(4,25),  THPIGZ(4,6), 

5E  RACT (4,6,25) 

COM  MON /DENS /TE , NO ALT, 1ALT( 10) , XUENS ( 4 ) , XFDE N 
COMMCn/GRIC/JNTGRID,XGRID(400) , VGRIDC400) ,PRCGRID(400) 

COMMON  /POINT /I  PCI  NT (4 ,20) , II PONT (4) , IDPONT (4  ) 

COMMON/ ALPHA/ NAME (4 ) .STATE (4,20) , 1ST  ATE (4 ,6) 

WRITE (2,70) 

70  FORMAT  (/6H1*****,*  THE  PROBLEM  CONDITIONS  ARE.*/) 

WRI1 E (2 , 7 1 )  Tt 

71  FORMAT(F10.2,*  =  E1ECTR0N  TEMPFHATURE  *) 

WRITE (2, 72)  X  EDEN 

72  FORMAT  (E10.2,*  =  ELECTRON  DENSITY  (CM-3)*) 

WRITE  (2, 62)  (C  IE  CC MKt)  .NAME  (ICE)  ,  ICE»  I  ,  NOSPEC  ) 

62  FORMAT  (F10.2,*  s  BEATY  S  IONIZATION  EQUATION  CONSTANT  FDR  *,A10) 
WRITE  (2,73) 

73  FORMAT  (/•  THE  HEAVY  BODY  NUMBEH  0EN5ITIE5  ARE.*) 

WRITE (2, 74)  ( NAMfc ( KK ) ,KK*1 .NOSPEC) 

74  FORMAT  (10(2X,A10) ) 

WRITE ( 2 ,75) (XDEnS(KK ) , K  F  c 1 .nCSPEC) 

75  FORMAT  (1UE12.3) 

NINT«INTGRID-1 

lalBGN 

D0  171  J»l,I 


PPDFA(J)«PROGRIO(J) 

371  CONTINUE 
IFLAG*0 
ICNFG«0 

CFPRQ»P  RDF A ( 1 ) 

CE5»(XGRID( 1+1 )-XGR10l I) 3/10,0 
ELESCsO.O 

OFPPOSU  )»0.2*DFPRO 
CFPROS(10)«DFPROS( 1 ) 

DFPROS(2)*0.4*OFPPC 

DFPR0S(9)«DFPH05<2) 

DFPR0S(3)«0,8*DFPRU 
0FPR05 ( 8  )«CFPRGS( 3 ) 

DFPPOS(4)«1.6*DFPRC 

DFPR0S(7)*CFER0S(4) 

CFPROS(5)s2.0*DFPRO 

DFPR0S(6)=CFPH0S(b) 

CO  15  J*  1 » NOSPEC 
KKN AXbNOFXC ( J  ) 

CO  16  K*  1 , K  NH  A  X 
PPGBE(J,K)r 0.0 

16  COMlhUt 
XK«AXeMONFS(  J3 
DO  17  K  =  1,KM<AX 
PRCBI(J,K)*0.0 

17  COM  INUE 
15  CONTINUE 

DO  14  J»l,10 
PRDFASC J)=0.0 
14  CONTINUE 
DO  2  0*1,10 

C  CLEAR  COLLISION  fPEQUENCY  ARRAYS 

CO  3  K  =  1  ,  NCSPFC 
KKMA  =  M0NFS(R) 

CO  4  KKsJ,KKRA 
CF IlK,KK)=0.0 

4  CONTINUE 
NKMX  =  N0FXC(><3 

CO  5  KK*1,KKHAX 
CFFCK , KKlsO.O 

5  CONTINUE 
3  CONTINUE 

C  INITIALIZE  POINTERS  FOR  EXCITATION  AND  IONIZATION  GRIDS 
CO  1000  K»l, NOSPEC 
KKNAXsNOFXC(N) 

IIRCNTCK)*! 

1 DPONT ( K ) a  1 
CO  1001  J  J  Js  1  , K  K  M  A  X 
IP01NT (N,0JJ)=1 
1001  CONTINUE 
1000  CONTINUE 

C  CALCULATE  THE  ENERGY  FOR  THE  BIN,  ALSO  THE  VELOCITY 

EsXGRIOU  +  l  )*0,S*DFS-F'10AT  ( J)*CES 
V»5,93E*7*SQRTIE) 

CALL  DEDTtF ,Tt,XELEN,DET,XL) 

ELCF  *-CET/CES 
SUMtXsO.O 

C  CALCULATE  THE  EXCITATION  COLLISION  FREQUENCIES 

DC  6  N  *  1 , NOSP  EC 
RKHAXsNOEXC (X  ) 

00  7  XK*1,KKRAX 


IF(E,LT,THRESH(K,KK))  GO  TO  7 
IFfJ.NE.l)  GO  TO  26 

25  LI«*1 

GO  TO  27 

26  LI*0 

27  CFE(K,KK)aOEX{K,KK,F,LI)*V*XDENS(K) 

SUHEX*SUNEX*CFE(K,KK) 

IFLACbJ 

7  CONTINUE 
6  CONTINUE 

CALCULATE  THE  IONIZATION  COLLISION  FREQUENCIES 
LOCATE  ENERGY  INDICIES  FOR  FRACTION 
CO  6  K«  1  , NOSPEC 
IF(E.LE.EICNGD(K, 1) )  GO  TO  8 
1 1 1 » 1 

10  IF(E.GT.FICNGD(K,lII).AND.F.L£.EIONGD(K,III+in  GO  TC  9 
IllBlIItl 

GO  TO  10 
9  UJCKJbIH 

8  CONTINUE 
SUNIONsO.O 

CO  11  K*1, NOSPEC 

oibqion(k,e  ,m 

RXNAX*MONFS(K) 

CC  12  KKcl'KKMX 
IKE.LT.THRlCZCK.KK))  GO  TC  12 
IIJKK*IIJ(K J 

CFIlK,KK)=GICfi(N.t,LI)*V*XCENS(K)*FRACnK,KKf  1IJKK) 
SUNION*SUMICN+CFI (K,KK) 

10NFG*1 
12  CONTINUE 

11  CONTINUE 
SCF*SUNI0N*SUMEX+EICF 
PROBIONaSUNION/SCF 
PR0BFX«SUMEX/5CF 
PPOBEL=ELCF/SCF 
PR00*DFPRGS(J)>PP0FA5(U) 

XNNX*PPOO*CES/SCF 
XNFE*PPOO/SCF 
FEL$(J)=V*PRCU/SCF 
IF(O.EG.IO)  GO  TC  107 
FR0FAS(J4  1)*ELCE  *(PRCO/SCn 

107  CONTINUE 

ACCUHILATE  PROBABILITIES 

ELESC  IS  THE  ENERGY  LCSS  TO  ELASTIC  SCATTERING  PEP  ELECTRCN 
FlE5C»ELESC«-(EtCF*DES*XNNX) 

CO  18  K»t, NOSPEC 
KKHAXbNOFXCCK) 

CC  19  KKsl.KKMAX 

PHCBE(K»KK)aPHOBFCK#KM'tCEE(K,KM*XNNX 

19  CONTINUE 
KKMAXbNIONFS(R) 

DO  20  RRsl.KKKAX 

PROBI  (X  ,KK)apROBl  (K  .MO+CFI  (K,HK)*XNNX 

20  CONTINUE 
18  CONTINUE 

LOCK  AFTER  DEGRADED  ELECTRONS  FRO*  EXCITATION 
XXBXNEF 

CO  21  Fal, NOSPEC 
XMAXaNOEXCm 
00  22  RXal , KKRAX 


Ifr'U. LI. THUtSHU, A*))  GO  TC  22 
FE«F-TRfiFSHlA,KR) 

A  A  A»  1 

24  IF(H.GT.XC.RlD(AAA>.ARr.Ft.LT.XGRIC(KAAfl))  GC  TC  2J 
KRK«  AR  A  1 

GO  TC  24 

23  PHDFA(RKK)BpRDFA<RKK)*AX*CFF(K,Ah)4lCtS/(AGRJD(RR|<4l)-XGRIC(KKII})) 
22  COMIM/t 
21  CUNTIMIF 

LOCK  AF7F.P  I'F  GKAl'Ft  ELECTRONS  FRO*  1CMZAT1CN 
CO  1  22  A«1,MSPFC 
XRMAX«MONFS(A) 
rc;  2»  a  A  « 1  ,  r  A  m  A  x 
lME.LT.THMCZlA.KR)  )  GO  TC  2H 
I1JAK*I1J(A) 

7F>P«XX*XI)FNS(A3»FRACT(R,AA,1IJAK)*V 
F SR AX»F -TMP IC21A ,KR ) 

AAhsl 

3i>  IF  (FSAAX.CF .XGR IT ( RRA ) . ANO.F S*AX ,LT. XCRIDCK AK+ 1  ) )  GO  TC  29 
RrRsRRK*  t 

GO  10  30 

29  FSPAXxO.MFbPAX 
RRG=  1 

<2  IF  l  F  S«AX  ,GF  .  XGR1CCRRC)  ,  ANO.FSMAX.LT  .  Xr.MCCRRCt  1)  )  GO  TC  31 

RKOxMvt*  1 
GO  TO  32 

31  IF  ( RRR ,GT.  I)  GO  TO  3  3 
F  HI  =2  .G*F  SP  AX 

I  F(0F  A(  1  JaF’Hl'F  At  IH2.0*U  Mp*oiF  IOMK,F  ,1  l,O,0,FUL,dJ*PY  )  *  IDES  /  I 

l XCR 10l 2 )-XGK  I  !H  m  ) 

GO  IC  2H 
U  ROKlxRRG-1 

00  34  A  R I  *  1  |RUR| 
l F ( RRC  .  F  U.  U  C.O  VO  3  9 

F’  K  0  F'  A  t  A  R  t  )  *  F‘  N  0  F  A  (  R  R  1  )9lFMF’4CIFI0N(R*F  ,  l  l,XGHIC(RRl  )#XGH1C(RK1f1)»C 
1 UMPY  )  •  (  DFS  /  (  XGKIC  (  R  K  1  1  3-XGMO  (  RR  1  )  )  ) 

3  4  COM  I  NO F 

IF  l  R  R(J.F  O.  1  )  GO  10  |  0  3 

3'kDFAIRRCJ  )  SF'WOF  A  (  KRO  ) -U  tPO’DlF  U M A  ,F  ,  l  I  ,  XGR lU( RAO  ),FSPAX,CU 
1  RR  Y  )  * l  T  F  S  / l XGF IOIrAOFI )  •  X  G  Ft  ID(RRC)  )  ) 

I  o  3  l  OM  I  Ml F 

F.SR’AX  =  2.0*F.S»*AX 

OF  PI sF  SP AX-XGA IlH  AAR  ) 

FHnFA(RRR)«F’KI>FA(RAR)-*lpMP*ClFlON(R,F,ll,0.0#UPPl 
I  ,  l  UM»  Y  )  *  (  OF  S/  l  XoP  l  (  (  R  HR  ♦  1  )-xGH  I  0  (  F  RR  )  3  ) 

IF  (RRu.FO.l  1  Gl  I V  I  04 
ULsFSRAX-XoRII  iRM.  ♦!) 

OF  PL*  0,9*  ES  "A  X 

PROF  A  l  RRG  )  =  PPOF  A  (  ARo  )  *T  F  *P*OI  MGA(R  ,  t  ,  l  1 , 01  l  ,  UPPl 
1  , 00  P  P  Y  I  *{CF  S/(XGNir(RRC*I  ) -Xf.h  I  0  I  RRli  )  )  ) 

104  COMIM'F 

FFU*Xs«Mi.| 

►  R  (.'  M  IS  S  R  R  0  ♦  1 
n  39  RR  1  *RRC»  r. ,  RRG«X 
I  F  (  RRG**N  .GF  .  RMJP  X  )  GO  TO  39 
ULl *F S* AX-XGH  1 1  IRRl  +  l) 

I  PFl*FSPAX”XOF>ll  (RRl  I 

FK0FAtRRI3  =  PRIFACRM)»TFMF*lIFll.'MR,F,lI,Ull,0PF'l  ,C 

1  U  P  P  Y  )  ♦IDFS/IXgKIC(RR1«1  )  •  X  G  K  1 1  I R  R  1  )  )  ) 

39  COAT  IM  F 
2  H  CGMlRUF 
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127  CONTINUE 
?  CONTINUE 
SUNsti.O 
CO  36  Jst ,  1  0 
SUN«SUN4FEIS( J> 

36  CONTINUE 

F EL ( I ) =  SUN  / 10. 

IF(I.EO.l)  GC  1G  37 

PRDFAt 1-1 )*PH)FA(I-1 )*tLCF*(PKnO/SCF)*(DFS/CXGPlC( I)-XGRIC(1»:  )  )  ) 
INN*I-J 

CG  38  IPs  1  ,  1 N  N 
IPVal-IP 
I N  V  *  PIN  NUNbtK 

CLEAN  collision  FREQUENCY  ARRAYS 
CC  JR  Ka  1  , NUSPEC 
CO  4'>  K Ns  1 ,6 
CFI(K,KK1*0.0 

40  CONTINUE 

K  Kw  AXaNOF  XC  (N  ) 

CO  41  KKal,KKPAX 
CEFlh.KNJxO.O 

41  CONTINUE 
38  CONTINUE 

CE=XGKI0( IRV+l J-XGRID(IRV) 

E=XfcRlC(IRV)+0.5*Gfc 

V*5,83E+7*SOPT(E) 

CALL  DECT  (  F  ,  1 1  ,  XEUE'N  ,  UE  T  *  XL  > 

E LCF=-CET/Ct 

CALCULATE  1HE  EXCITATION  COLLISION  FREQUENCIES 
SUNEXsO.O 
CC  42  Kaj, NUSPEC 
KisfcAXsNOEXC  (X  ) 

Cl)  43  hKsl,Mt*AX 
1ME.LT.THPESm(|*,nNJ)  GC  TC  41 
CFFlK,NK)aCtX(K,KN,E,lt)*V*XDENS(K) 

SUMFXaSUMEX+CFECK.KN) 

43  CONTINUE 

42  CONTINUE 

CALCULATE  THE  10NIZAIIGN  COLLISION  FREQUENCIES 
LOCATE  THE  ENERGY  IMUC1ES  FOR  EKACTICN 
ro  44  K=1 , NLS  P  EC 
IE ( t.LE .FICNCL(K , 1 ) )  GC  TO  44 
I  11  =  1 

45  IE  (E  .GT.EICNC  0(N,  III)  .  AM)  .  F  .  LE  .  F.  IGNGC  (  K  ,  1 1  !♦  1  ))  GO  TC  46 
IIIaIiI+i 

GC  TO  45 

46  IIJlKjaltl 

44  CONTINUE 
SUNIONaO.O 

CO  47  K  a l , w OS PEC 
NK*AX  =  MONE  S(  N  ) 

CC  48  N Ka 1 ,KM  AX 
IF(E.LT,THR1CZ(N,nX  )  )  GO  TC  48 
IIJKKaflj(K) 

CFI(K,XF)aUlCN(F ,E , L I) » V ♦ XU t NS ( M *F * ACT ( K , KN , 1 1 JUK  ) 
SUNIONsSUNlONtCFl CK.KK) 

48  CONTINUE 

47  CONTINUE 
SCFaELCF^SL REX^SI  HON 
PKCPIONaSOUCN/SCF 
PROBF  X  =  SUHE X/SCE 


<>?) 


PFtOb£L»ELCF/SCF 
GRQ*PRDF A( IRV  } 

XNNX*GR0*Dfc/5CF 
FEL(IRV)*V*GRQ/SCF 
XNEE*GRO/SCF 
IFlIRV.EQ.l)  GO  TO  38 

PRDFA(IRV-1)*PRDFA(IRV-| )+ELCF* (GRC/SCF)  * ( CE/ ( XGR I D ( IP V ) -XGRIC C 

1IRV-1))) 

C  ACCURJLATE  PROBABILITIES 

FLESC=£LESC4(ELCF*DE*XNNX) 

CU  49  K=1,N05PEC 
KKRAXaNOEXC(K) 

CO  50  KK  si ,  K  K  R  A  X 

FRCBECK  ,KK)=PKOBElK,KK)4CFE(K,MO*XNNX 

50  CONTINUE 
KKNAX  =  MONFS(K) 

DO  51  KKsl.KKwAX 

PROBI  (K,KK)sPROBl(K,KK)+CFnK,KlO*XNNX 

51  CONTINUE 
49  CONTINUE 

C  LOOK  A 1  TEH  DEGHADED  ELECTRONS  F'RC*  EXCITATION 

XXsXNEE 

00  52  Ks  1  »  NCjSPFC 
KK*AX=NOfXC(K) 

DC  53  KKrl.KKMAX 

IF (E.LT .THRESHCK ,KK) )  GO  TC  53 

fcE=F>THRFSh(K,KK) 

KKK  =  1 

55  IF  (FF:.GF.XGR1DCKKK),ANP,EK,LT,XGRIC(KKK  +  t))  GC  TO  54 
KKKsKKKM 

GO  TO  55 

54  PHDFA(KKK)»PRDF A ( K K K ) ♦ XX *CFE ( K , KK ) * ( CE/ ( XGR IC ( K KtU  1 ) -XGR IC ( KKK ) ) ) 

5 3  CONTINUE 

52  CONTINUE 

C  LOCK  AFTER  DEGRADED  ELECTRONS  FRO*  IONIZATION 

CO  56  k  =  i,nospec; 

KkMAXsNIGNFS(K) 
no  57  KKal.KKMA* 

IME.LT.THR1CZCK.KK1  1  GO  TC  57 
ESMAX=F-THRIGZ(K,KK ) 

IlJKKsIIJ(K) 

TERPsXX*XOENS(K)*V*KHACT(Kf KK, UJKK  ) 

KKK=1 

59  IF CFSHAX .GF.XORID (KKK ) . AND.ESMAX.LT.XGPICC KKK+1) )  GO  TC  56 
KKK=KKK-f  1 

GO  TC  59 

56  ESNAX  =  «I.5*ESRAX 
KKO  =  l 

60  IF ( ES*AX .Gfc .XGR1C (KKO) . AND.ESMAX.L7 . XGPIC (KKC+l ) )  GO  TC  61 
KKOsKKO-M 

GO  TC  60 

61  IF (KKK.GT.il  GO  TC  6b 
tUL=2.0*FSR AX 

PR0FA( 1 JsPROF A(l ) ♦ 2 . 0 *T E*P * 0 1 F I ON ( K , t , l 1 , 0 , 0 , EUL , DUMP Y )  *(DE/(XG 

1PIOC  2J-XGR  I  D  C 1 ) ) ) 

GO  TO  57 
66  KUMSKKU-1 

CU  67  KK1  =  1,K0M 
IFCKKO  ,  FO  •  1 )  GO  TO  6  7 

PROFACKM )»PPDF  ACkKI ) 4TE"P»D IF ION ( A . E , l I , XC R I C ( K K 1 ) , XGR I C ( KK 1 ♦ 1 ) , C 
1UMMT)  41CF /IXGRIDCKKI+l J-XGRIDCKFl ) ) ) 
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67  CONTINUE 

IF(KKQ,E0.1)  GO  TO  105 

FRDF  A  (  KKO  )«PRUFA(KKO  ) fTEMPFOIFION ( K , E , L l , XGRI D (KKQ  ),ESMAX,DU 
1  MMY  )  *(Ct  /  (  XGRIDUKQtl  )-XGRID(KKC  )  )  ) 

105  CUNTINUE 
ESMAX«2.0*ES»*AX 
UPPL*ESMAX»XGR1D(M(K  ) 

PROFA(KKK)*PRDFA(KKlO*TEMP*OIFION(K,E  ,H,0.0,UPPL 
1 .DUMMY )*(DE  /(XGRIOtKKK'Fl  )-XGRID(KKK)  )  ) 

1F(KKQ,EQ,1)  GO  TO  106 

ULL«ESMAX-XGHID(KN0+1  ) 

UPPL«0.5*ESMAX 

PRDFA(KKQ)aPRDFACKKO)+TEMP*DIFION{K,E ,LI,ULL,UPPL 
1  .DUMMY  )  *(DE  /(XGRIDCKKC+1)*XGRID(KK0))) 

106  CONTINUE 
KKQ  MXaKKK  *•  1 
RKQMNsKKG* 1 

DC  68  KK1«KKQMN,KKCMX 
IF(KKQMN.GE.KKQMX)  GO  10  68 
ULL  =  ESMAX-XGPID(KM*1 ) 

UPPL=ESMAX-XGRID(lvKl ) 

FHDFACKK1  )  apR OF A (  M<  1 ) M EMP* 0  IK  ION ( K , E , L I , UL L , UPPL  ,C 

1UMMY)  *(Dt/(XGRlD(KKl*1  ) -XGRID ( KK 1 )  )  ) 

68  CONTINUE 
67  CONTINUE 
56  CONTINUE 
38  CONTINUE 
37  CONTINUE 

Es(XGRID(l+l)+XGRID(l))*OaS 

Va5.93EF7*SUHTCfc) 

WHITE (2, 76) 

76  FORMAT  C/6H  ♦  ♦*♦♦,*  SUMMARY  OF  EXCITATIONS  PRODUCED.*) 

IF (IFLAG.EC.O)  GO  TC  7 7 

SUM  1 »0 , 0 

CO  78  K  =  1 »  NOSREC 
MRITE  (2,75) 

WHITE(2,74)  NAMt(K) 

WRITE(2,75) 

KKMAXsNOFXC  (  F  ) 

MRITE(2,74) (STATE (K,KK ) ,KN=1 ,KNMAX) 

WR1TE(2,75)  (PRUPt(K,KK),NK=l,KKMAX) 

SUMsO.O 

DO  86  KKeV ,  K  K  M  A  X 
5UM*SUM*PRCbF ( K , KK  ) 

86  CONTINUE 

WR1TE(2,87)  SU«,NAMt(K) 

SUMlSSUMl+SUM 

78  CONTINUE 
WRITE(2,88)  SUM1 
GO  TO  79 

77  «RITE(2,80) 

80  F0RMAT(//,69H  ti'F'RGY  IS  HELD*  ALL  EXCITATION  THRESHOLDS,  LOSS  TO 
1  ELECTRONS  OMY  ,//) 

79  CONTINUE 
WR1TE(2,81 ) 

81  FORMAT  C/6H  *♦*•*,♦  SUMMARY  OF  IONIZATIONS  PRODUCED.*) 
IFdCNFG.EG.O)  GO  TC  8? 

SUMlso.O 

CO  83  Raj , NOS  PEC 
mRITE(2,75) 
wRITF(2,74)  NAME(k) 
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hRITE(2,75) 

KKMAXaNlONFS(K) 

*HITE(2»74)  (ISTATE(K,KK) ,KKal,KKMAA) 

NRIIE(2,75)  (PR08I(K,KK) ,KK*1 ,KHMaX) 

SUMsO , 0 

00  89  KK*1,KKMAX 
SUMaSUM*PRCBI(K,KK) 

89  CONTINUE 

*RI1E(2,90)  SUM,NAMF(IO 
SU»1*5UN|+SUN 

83  CONTINUE 
SPACEsO.O 

AC202aPR0BI( 1 .1 )*PRGBI(l,3)+PROBI(l,5) 

AC402*PR0BI(1 ,2)  +  PR0BI(  1 ,4) 

ACN2»PR0BI(2, l ) +PROB I ( 2 , 2 ) +PROB I (2,3)*PR0BI (2,4) 

ACNP*PROBI(2,6) 

b RITE (6,926)  ACN2  ,  AC202 , AC4C2 , SPACE , PRCB I ( 3 , 1 ) , PROB I ( 3 , 2  )  , 

1FR0B1 (3,3) ,ACNP 

926  FORMAT  (7X.6E10.3) 

ELESC»ELESC+(5UN1*T0TF  L ) * ( XGR 1 0 C 2 ) -XGR ID ( 1) ) 
bRITE (  7 , 927  )  ELESC 

927  FORMAT  (8E10.4) 

*RIIE(2,91)  SUM  1 
WH1TE(2,110)  ELtSC 

110  FORMAT  (//F10.2,*  =  ENERGY  LOST  TO  THERMAL  ELECTRONS  (EV  CM-3* 
t  ,*  SEC-1)*) 

87  FORMAT  (/  tlo.3,*  *  TCTAL  EXCITATIONS  OF  *,A10) 

88  FORmATC//  E10.3,*  s  TOTAL  EXCITATIONS*) 

9C  FORMAT  (/  t 1 0 , 3 , *  *  TCTAL  IONIZATIONS  OF  *,AlO) 

91  FORMAT (//  E1U.3,*  a  TCTAL  IONIZATIONS*) 

GO  TO  84 
82  bRITE(2,85) 

85  FORMAT!//, 81H  ENERGY  IS  BELCm  ALL  1LN1ZAIICN  THRESHOLDS,  LOSS  TO 
1  EXCITATIONS  (  IF  POSSIBLE  )//,17h  AND  ELECTRONS  ) 

84  CONTINUE 
*HITE(2,93) 

93  FORMAT!//*  EQUILIBRIUM  FLUX  DISTRIBUTION  (CM3  SEC  STER  EV)-1*) 

CC  94  FK*I,1 

FEE(KK)a(XGKI0(l\K)*XGHlD(KhEl))*0.5 

94  CONTINUE 
*RIIE(2,9b) 

96  FORMAT  (/b(9X,lHF,9X,IHF)/) 

653  FORMAT (8E10, 3) 

DO  945  KK=1,I 
FEL(KK)aFF.L  (MO /12,56b 
945  CONTINUE 

*R1TE(2,97)  ((EEE(KK),FEL(KN)),KK  =  1  ,J) 

WRIT  t  (2,2234  )  1AVGF , IAVGF 
IF  (IAVGF. IT, 2)  GO  TO  6314 
00  2227  KK=J,1 
PROFA(M  )  =  0.0 

2227  CONTINUE 
lAveaal AVGF/2+1 

OC  2228  KEal, IAVGF 

PROF  Ad  JaFRDF  A(1  )*F'FL(KK) 

2228  CONT 1 NUF 

PROF  A ( 1 )=PHDFA( 1 )/F  LOAT(IAvGF) 

00  2229  EKa2, 1AVBH 
PHOFA(KK)aPROFAd) 

2229  CONTINUE 

I AVCCal-IAVBrt 
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1 A VAA*1 A VBB4 1 

DO  2230  KK3IAVAA,1AVCC 

ULTsKK-IAVGF/2 

IVLTBKK4IAVCF/2 

DO  2231  1 J AVsILLT  ,  I  V  LT 

PRDFA(KK)*PPDFA(KK)4FEL(IJAV) 

2231  CONTINUE 

PRDFA(KK)»PR0FA(KK)/FL0AT(1AVGF) 

2220  CONTINUF 

IAVCCSI AVCC41 

DO  2232  KN*IAVCC,  I 

PRCFA(1AVCC)*PRDF»(1AVCC)4FEL(KK) 

2232  CONTINUE 

PRDFA(1AVCC)=PKCF A ( 1  A VCC ) /FLU AT (1 A VGF  ) 

IAVCC=IAVCC+1 
00  2233  KK=IAVCC,1 
PRDFA(KK)  =  PRDF A(IAVCC4l  ) 

2233  CONTINUE 

2234  FORMAT  (//I3,*  (LUX  PEP  STERAD1 AN  AVERAGED  OVER*, 13,*  INTERVALS*) 
*R1TE(2,96) 

RRME(2,97)((F.EE(KK),PRDFA{KK)),KK*1,I) 

6314  »P 1TE ( 5 , 946 )  1 

946  F0RMAH15) 

*RITE(5,947)  C FEL ( K K ) , KK* 1 , 1 ) 

•  RITE ( 5 , 94  7  )  (EEE‘(KK),MS1,1) 

947  FORMAT! 6E 12.4) 

97  FURMAT(6(F1U,3,E10,3)  ) 

IF  HAVGF.LT.2)  GO  TO  6316 
LC  900  RK*  1 , I 
FEL(KK)*12,568*FEL(KK) 

900  CONTINUE 

RRITE(2,30P) 

300  FORMAT  (//  *E0U1LIBRUIM  FLUX,  NOT  PFR  STERALI AN* ) 

*RITE(2,97) l (EFE(KK) ,FEL(KK)  )  ,K6«1  ,1) 

6315  RETURN 
END 

SUBROUTINE  REDXSN 

C  THIS  ROUTINE  READS  THE  CRUSS  SECTION  DATA 

COMMON/XSXS/  NOSPtC,  NCEXC(4),  nGRIC(4,20), 

1EGR1D(4,20,25),SPLCO(4,20,26),SPLC1  ( 4 , 20 , 25 ) ,SPLC2 ( 4 , 20 , 25  ) , 

25PLC  3(4,20,25)  ,  THRESH  (  4 , 20  )  ,  MONFS  (  4  )  ,  MCNC.L  (4  )  ,  EICNGD  (  4 , 25  )  , 
3CIFCCNC4) ,T 1 X  S  0  t  4 , 2  5 ) ,TIXS1(4,25),TIXS2(4,25),TIXS3(4,25),CIXS0(4, 
425) ,D1XS1 (4,25) ,DIXS2(4,25) , 0 1 XS 3 ( 4 , 25  )  ,  7RR1CZ(4,6), 

5FRACT (4,6,25) 

COMMON/ALHHA/NAHF (4) , ST  ATE (4 ,20) , 1ST AT F( 4, 6) 

DIMENSION  1FCRMC8) 

REACH  ,20)  ( 1  FORM  (I),  1  =  1 ,8) 

20  FORMAT(IOAB) 

REACH, 1)  NOSPEC 
l  FORMAT ( 2014  ) 

C  NOSPEC  IS  T  .'t  NUMBER  OF  SPECIES  PRESENT,  LIMIT  4 

CO  2  1*1 ,NCSPFC 
RE AD ( 1 , 3  )  NAME(1),N0EXC(I) 

3  F0RMAT(A10, lb,F6.2) 

NAME(I)  IS  THE  NAME  OF  THE  SPECIE  I 

NOEXC(I)  IS  THE  NUMBER  OF  EXCITAUCN  CROSS  5ECTI0N5  FCR  SFEC1  1 
NsNDEXC ( I ) 

CO  4  0*1,  N 

REACH, 3)  STATE!  i,J),NGPID(l,J),THPESH(  1,0) 

STATE  (1,0)  IS  THE  STATE  CE5IGNAT1CN 
NGPID  (1,0)  IS  THE  NUMBER  OF  POINTS  IN  THE  GRID 
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C  THPE5H(I,J)  IS  1H£  THRESHOLD  ENERGY 

N1 =NGRID ( 1,0) 

READ (  1 ,  l )  l.FCKH 
IF (LFOPR ,E0 , 0  )  GO  TO  21 
READ  (I  .irOPMHEGRIDCI  ,0,K)  ,K*l,Nl ) 

GO  TO  22 

21  PEAD(1,5)(EGP1D(I,J,K),M1,M) 

22  R1«M«*1 

REAU(l,6HSPLC0(I,0,K),K*l,Nl) 

READ  ( 1, 6  )  (SPICK  1,0,  K),K«1,M) 

READ ( 1 ,6) (5PLC2(I ,  J  ,  K  )  ,  K*  1  ,N1  ) 

REAOC1 ,6) <SPLC3(I , J,K)  ,K«1  ,N1) 

READ (1,1)  KYY 
lE(KYY.EQ.O)  GO  TO  4 
READ (1,6)  CORCT 
CO  9  M1,M 

SPLCO(I,J,K)*SPLCO(I,J,K)*CORCT 
SPLCK1,  J,K)*SPLCKI,J,K)*CORCT 
SPLC2CI, J,K)»SPLC2(I, J,K)*CORCT 
SPLC3 ( I , J , K)«SPLC3( 1 , J , K ) *CQRCT 
9  CONTINUE 

5  F0RNAT(10C1X,F6,2)  ) 

6  F0RM1(7(1X,E  10.3)) 

EGR10  IS  THE  ENERGY  GRID  FOR  THE  SPLINE  FITS 
SPLC0,SPLC1,SPLC2,SPLC3  ARE  THE  SPLINE  COEFFICIENTS 
4  CONTINUE 

RE AC ( 1 , 1 )  NIONFSf I) .NlCwGDd) 

REA0(1,S)  DIFCCN(I) 

MONFS(I)  IS  ThE  NUMBER  OF  ION  FINAL  STATES, LIMIT  6  FOR  EACH  SPECIE' 
NIGNGO(I)  IS  THE  NUHBER  Of  GRID  POINTS  FOR  IONIZATION  CROSS 
SECTION  DATA 

DIFCON(I)  IS  THE  CONSTANT  IN  BEATY ) 5  IONIZATION  EQUATION 
N2*NI0NFS(  l ) 

N3«M0NGD(  I  ) 

RE AD (  1 ,  1 )  LEFORm 
IF  (lEFURR.fC.O)  GO  TO  23 
READ  ( 1 ,  IFORn  )  (EIONGDd  ,J),J*1  ,N3) 

GO  TO  24 

2 J  READ ( 1 ,6) (EICNGO(l.O) ,0=1 ,N 3) 

24  N3=n3-1 

REAu(1,6)(T1XS0(I,J),J=1,n3) 

REA0(l,b)(IlXSKI,0),J»l,N3) 

Ri A0( 1 ,6) (TIXS2( 1,0) , J*l,N3) 

REA0(1,6)(T1XS3(1,J),J«1,N3) 

REA0(l,6)(CIXS0(I,0),J«i,N3) 

R  t  A  D ( 1 ,6) (C1XS1 (I ,J) *  J  =  1 , N  3 ) 

READ(l,b)(DIXS2(I,J),J=l,N3) 

RE AO (1,6)  (C  1X53(1,0)  ,J=1  ,M) 

REAC(l,l)  KYY 
IF (KYY.EO.O)  GO  TO  10 
REAO(l.b)  CORCT 
CO  11  0=1 ,N3 

TIXSO(l ,J)*TIXS0(I,0)*CURCT 
1IXSK1,J)*TIXS1(1,0) ♦CORCT 
TIXS2(  I,0)«TIXS2( I,0)»C0RCT 
T1XS3(I,O)s1IXS3(1,O)*C0RCT 
D I  X  50 (  I,0)*DIXSu(l,0)*CORCT 
C 1XS1 (  1,0)=DIXS1(1,0)*C0RCT 
LlXS2(l,0)=DIXS2(l,0)*C0RCT 
D  I  X  S  3 ( 1 ,0)  =  DIXS3( 1 ,0 )*CORCT 
11  CUMINUt 
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10  CONTINUE 

E10NGD  IS  THE  IONIZATION  CROSS  SECTION  ENERGY  GR10 
TIXSO  TO  TJXS3  ARE  THE  TOTAL  IONIZATION  CSOSS  SECTICN  SPLINE 
COEFFICIENTS 

OIXSO  TO  C1XS3  ARE  THE  DIFFERENTIAL  CROSS  SECTICN  SPLINE  COEFFICIENTS 
FOR  THE  NUMERATOR  IN  BEATY  CS  EQUATION 
DO  7  J«1,N2 

READ ( 1  ,  B  )  ISTATECl, J)  ,THR10Zd  ,  J) 

B  F0RMAT(AB,F6,2) 

1ST  ATE  IS  THE  FINAL  ION  STATE  DESIGNATION 
THRIOZ  IS  THE  THRESHOLD 
READ(1,5)CFRACT(I,J,K),K*1,N3) 

7  CONTINUE 
2  CONTINUE 

FRACT  ARE  THE  BRANCHING  RATIOS 
RETURN 
END 

FUNCTION  QIONCI ,E,KEY) 

COMHON/XSXS/  NOSPEC,  NCEXC (  4 )  ,  NGRID (4,20), 

1EGRID(4,20,25) , SPLCO ( 4 , 20 , 25 ) , SPLC1 C 4 , 20 , 25 ) , SPLC2 ( 4 , 20 , 25 ) , 
2SPLC3(4,20,25) , THRESH < 4 , 20 ) ,NI0NFS(4) ,  MONGD  (  4  )  ,  E 10NC0  (  4 , 25  ) , 
3D1FCCN(4),TIXS0(4,25),T1XS1(4,25),T1XS2(4,25),TIXS3(4,25),CIXS0(4, 

4  25)  ,D1XS1(4,25),C1XS2(4,2S) ,CIXS3{4,25),  THRIOZ (4, 6)  , 

5FR ACT (4,6,25) 

COMMON  /POINT/1  PC  I NT (4, 20) , I1P0NT(4) , IDPONT (4 ) 

COMMON /ALPHA/ NAME  14) .STATE (4,20), 1ST ATE ( 4,6) 

This  FUNCTION  CALCULATES  THE  TOTAL  ELECTRON  INPACT  IONIZATION 
CROSS  SECTION  FOR  THE  SPECIE  I  AT  ENERGY  E.  KEY  AND  1IFCNT  ARE 
USED  TO  SAVE  1IME  IN  THE  ENERGY  GRC  SEARCH,  IF  KtY»»l,  ME  START 
AT  THF  BEGINNING  OF  THE  ARRAY  IN  THE  SEARCH  AND  SET  IIPCnT  BEFORE 
EXIT.  IF  KEY  a  0,  nF.  START  AT  IIPCNT  AND  RESET  IIPOM  IF  NECESSARY. 
IE  KEY  IS  POSITIVE,  »E  START  AT  KEY  AND  RESET  IIPCNT  IF  NECESSARY 
120  FORMAT  (4H  I  3 , 1 4  ,  t>H  KEY  a, 14) 

110  F0RMATI4H  E  a,  E9.2.14H  EI0nGD(I,1)  a,  E9.2) 

1F(E.GF.EICNGD(1,1))  GO  TO  2 

1  QIONaO.U 
RETURN 

2  NnsnIONGC(I) 

1F(E.GT.FICNC0(I ,nn)  )  CO  TO  1 
JOO  FORMA T(5H  NN  =,I4,4h  I  =,I4,4H  E  a,E9,2, 

1 6 H  KEY  a, 14) 

3  CONTINUE 
lF(KFY)  4,5,6 

4  I  Pa  1 

9  IE (E.GE.FICNGUU , IP) , AND.t.LT.tlONGDt I  ,  1P+1  )  )  GO  TO  7 

8  lP*IP4l 
GO  TO  9 

7  1 1  PONT ( l )  a I P 

10  XaE-EIONGD(l.IP) 

YsElCNGO(I,IP+l)-i 

010NSX*CTlXS3d,IP)*TlXSld,IP)*X**2)*Y*(lIXS2d,IP)4TlXSU(l,IP)*Y 

1**2) 

RETURN 

5  IPslIPCNT(I) 

12  CONTINUE 

1F(E.GE.FICNGD(I,1P).AND.E.LT,FI0NG0(I,1P*1))  GO  TO  10 

11  IF(E.GE.FICNCU(1,1F))  GO  TO  8 

13  IP«IP-1 

1F(E.  GE.ElCNGOd,  IP),  AND.  t ,L1 . E 10NGD ( 1 ,  IP*  1 )  )  GO  TO  7 
GO  TO  13 

6  IP«KEY 
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IIPONT(I)»IP 
GO  TO  12 
END 

FUNCTION  OIF  ION U ,E,KEY,E0,E1 ,F) 

C0MM0N/X5XS/  NOSPEC,  nCEXC(4),  NGRID(4,20), 
lEGRID<4f20,25),SPLC0(4,20,25),SPLCl(4,2O,2S),SPLC2(4,20,25), 
2SPLC3(4,20,25)  .THRESH (4,20) , MONFSI4), NI0NCD(4), EIONGDU, 25), 
3ClFCON(4),TIXSO(4,2b),TlX51(4,25),TlXS2(4,2S),TIXS3(4,25),CIXS0(4, 
425),DIXSl(4,25),CIXS2(4,2b),0IXS3(4,25),  THRIQZ(4,6) , 

5FR ACT (4,6,25) 

COMMON  /P01NT/T POINT (4,20) , IIP0NT(4) ,  1 DPONT (4) 

COMMON /ALPHA/NAME (4), ST ATE (4, 20 ) , 1ST  ATE (4, 6) 

DIMENSION  F(6) 

THIS  FUNCTION  CALCULATES  THE  ELECTRON  IMPACT  IONIZATION  DIFFERENTIAL 
CROSS  SECTION  FOR  THE  SPECIF.  I  AT  ENERGY  E,  KEY  AND  IIFCNT  ARE 
USED  TC  SAVE  TIME  IN  THE  ENERGY  CRD  SEARCH,  IF  KEYs-1,  ME  START 
AT  THE  BEGINNING  OF  THE  ARRAY  IN  THE  SEARCH  A NO  SET  ICPGNT  BEFORE 
EXIT.  IF  KEY  s  0,  mE  START  AT  IDPCNT  AND  RESET  IDPONT  IF  NECESSARY, 
IF  KEY  IS  POSITIVE,  ME  START  AT  KEY  AND  RESET  IDPONT  IF  NECESSARY 

me  integrate  beatyjs  eolation  from  eo  tc  ei  using  a  cubic  fit 
mH  1  ch  IMPLIES  small  step  sizes.  The  fractional  yeald  to 

DIFFERENT  ICN  STATES  ARC  RETURNED  IN  THE  ARRAY  F. 

IF (E,GE.EICNGD(1,1))  GC  TU  2 

1  GION»O.U 
RETURN 

2  NNCNlONGD(I) 

IF(E,GT,FICNGU(I,NN))  GO  TC  I 

3  CONTINUE 
IF(KET)  4,5,6 

4  IP«1 

9  lF(E.GF,FXCNGD(l,IP),AND.E.LT,EIONGD(l,IFtl))  GO  TO  7 
8  IP  =  IP*1 

GO  TO  9 

7  IOPONT(I)*IP 

10  XsE-FIONGD(l.IP) 

Y=E10nGD(I , IP+1 )-t 

COhO*X*(DlXS3U,IP)MDlXSl(l,IP)MX**2)*Y*(ClXS2(I,IPHOIXSQ(I,lP)*Y 
1  4*2) 

MI  =  NI0NFS( 1  ) 

CO  14  J=1,ND 
F  (J)*FRACT(l,wl,IP) 

14  CONTINUE 

Em  1 0*0 ,  5*  ( E 1  +  EO) 

Di«0,S*(F.  I“t0) 

FXXsl.O/d.O+CEO/OIFCUNd)  )  *  ♦  2 , 1  )  4  4 , 0  /  (  1 , 04  (  E  M  I D/ D 1 FCON  ( 1  )  )  **2 , 1  ) 
141.0/ ( 1.04(8 l/OIFCON(l) )**2.1  ) 

C1F10N«C0NC*E  XX*Cfc/3,0 
RETURN 

5  IP*IDPOi*T(I) 

12  CONTINUE 

IF (E.GE.FlCFGDd ,IP) , AND . t . LT , EIONGC ( 1 , IPtl  )  )  GO  TC  10 

11  lF(F,GE',FlCNGD(i  ,  IP)  )  GO  TO  8 

13  IP*1P-1 

IF  (E.GE.E1CNGO(1,1P),AND.F:,LT,EIONGO(I,IP41))  GO  TC  7 
GC  TO  13 

6  IP«KEY 
IOPONT( 1)=JP 
GO  TO  12 
FNC 

FUNCTION  QFX(I,J,t,KFY) 

COMMON/XSXS/  NOSPEC,  N0EXC(4),  NGRID(4,20), 
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1EGR1D(4,20,25),5PLCO<4,20,25),SPLC1(4,20,25),SPLC2(4,20,25), 

28PLC  3 ! 4 , 20  <  25 ) , THRESH ( 4 , 20 ) , N IONFS ! 4 ) , N10NGD ( 4 ) , E1CNGD! 4 , 25 ) , 

30IFCQN ! 4) ,T1XS0 ( 4 , 25 ) , T IXS1 ( 4 , 25 ) , T1XS2 ! 4 , 25 ) , TIXS3 ! 4 , 25 ) , DIXSO (4 , 

425) ,0IXS1 (4 , 25) , C1XS2 ( 4 , 25) , DIXS3 (4 , 25 )  ,  THRI0I(4,6), 

SFHACT (4,6,25) 

COMMON  /POINT/IPC1NT(4,20),IIPONT!4) , IDPONT (4) 

C0MM0N/ALPHA/NAME(4) , ST ATE (4, 20) , IS) ATE (4, 6) 

THIS  FUNCTION  calculates  the  electron  impact  excitation  cross 
SECTION  FOR  THE  PROCESS  J  IN  SPECIE  I  AT  THE  ENERGY  E,  KEY 
AND  IPOINT  ARE  USED  TC  SAVE  TIME  IN  THE  ENERGY  GRID  SEARCH)  IF 
KFY  *-l,  HE  START  AT  1  HE  BEGINNING  OF  THE  ARRAY  IN  THE  SEARCH  AND 

SET  IPOINT  BEFORE  EXIT.  IF  KEY  aO,NE  START  AT  IPOINT  AND  RESET  IPOINT 
IF  NECESSARY!  IF  KEY  IS  POSITIVE,  HE  START  AT  KEY  AND  RESET  IPOINT 
IF  NECESSARY) 

IFtE.GT.EGRIDC I , J, 1) )  GO  TC  2 

1  0EX*0,0 
RETURN 

2  J  1  sNGR ID ( I , J ) 

IF!E.GE.EGPID!I,J,Jl))  GC  TO  I 

3  CONTINUE 
IF(KEY)  4,5,6 

4  IP»1 

9  1F(E.GE.EGPIDU,0,IP),AND.E,LT,EGRID{I,J,IFA1 ))  GO  TC  7 
8  IP*1P*1 
GC  TO  9 

7  IP01NT(I,J)*1P 

10  X»E-EGRID(I,J,TP) 

YsEGHIDII , J , 1  P*  1 )"E 

OFXaX*!SPLCJ(I,J. lF)tSPLCl( I, J , 1 P ) ♦* * *2 ) ♦ Y ♦ l SF LC2 ! I , J , IP)*5PLC0(I, 
10,IP)*Y**2) 

RETURN 

5  I  pa  1 PG I  NT ( I «  J  ) 

12  CONTINUE 

IFCE.GE.FGRIOCI ,  J,  IP)  .  AND  ,  E  .  LT  ,  EGR IDU  ,  J  ,  IF*  1  )  )  GC  TO  10 

11  IF  (E.GE.FGR1CU  .  J,  IP)  )  GO  TO  8 

13  IPsJP-1 
JF(IP.FG.O)  GO  1C  1 

IF(i.GE.EGRlC(I,J,lP).AND,E.LT.fcCRID(I,J,lF*l))  GC  TO  7 
GC  10  13 
t  IpaNEY 

IPCINTt 1 , J)*IP 

GC  TO  12 

END 

SUBROUTINE  SETGHIO 

SETGR1D  SFTS  CP  THF  GRID  FOR  INTEGRATION 
COMMON/GR 10/ INTGRIO.XGR 10(400) ,VGH 10(400 ),FHCGRID (400) 

DIMENSION  E0(10) ,EF(10) ,NPTS( 10) ,DFE( 10) 

RE AO ( 1 , 1 )  ITYPE 

ITYPE  SPECIFIES  THE  TYPE  Cf  GRID  DATA  .  IF  ITYPE  ai,  READS 
INTGRID  PCIMS  FOR  XGRIC)  IF  ITYPE  a2,  READS  NRANGE  SETS  OF 
INITIAL  VALUE SCFInAL  VALUES!  AND  NUMBER  OF  PClNTS 
IF  ITYPE  *3,  READS  INITIAL  VALUE  AND  NRANGE  SETS  CF  INTERVALS 
ANO  NUMBER  CF  PClNTS 

JTYPE«4,  XGR1U  IS  SET  IN  PRODUCTION  SUBROUTINE 

1  FORMAT (14) 

CO  TO  (2, 3, 4, 341), ITYPE 
341  RETURN 

2  RE AC ( 1 , 1 )  INTGRID 

READ  (1 ,5)  IXGMlDd  ),I»1 ,  INTGRID) 

5  FORMATUO!  IX.F'6.2)) 

CC  12  Kal ,  IMGRIO 
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r*  r>  n  n  r> 


r 


VGRID(K)«5,93i>7*5QRT(XGRlD(K) ) 

12  CONTINUE 
RETURN 

3  REAU(l,n  NRANGE 
CO  6  1*1, NRANGE 
READ (  1,7)  EOm.NPTStn.EMI) 

7  FORMATUX,F6.2,2X,  14,2X,F6.2) 

6  CONTINUE 

INDEX*! 

DO  8  1*1, NRANGE 
XGRID(1NDEX)*EU(I) 

Dt*(EF(l)-EO(I))/FLOAT(NPTS(I)-l) 

N  M  AX*NPTS  C I  ) 

CO  9  J*2,NNAX 
IND£X*lNCEX*l 

XGR ID ( INDEX )*XGRID( INDEX- ll+OE 
9  CONTINUE 

8  CONTINUE 
INTGRIDsTNCEX 
CO  13  K*1,IMGRID 
VGR I D ( F )=5,93E+7*S0RT (XGRID(K)  ) 

13  continue 

PEI  URN 

4  READ (1,5)  FUG 
REA 0(1,1)  NRANGE 
CG  10  1*1, NRANGE 
REAU(1,7)  CEFCI)  ,NP'IS(  l) 

10  CONTINUE 
iNDEXsl 
XGRID( 1 )=ECC 
CO  11  1*1, NRANGE 
NMAX  =  NPTS( I ) 

DO  11  J*1,N«AX 
INDEXsINDEX+1 

XGRID(INCEX)=XGRID(INOEX-l)iDEE(I) 

11  CONTINUE 
1  NTGM 10*1 NCEX 
DO  14  Ks 1 , I NTGR 1 D 
VGR10(K)*b.93t ♦74SORT(XGRIC(K) ) 

14  CONTINUE 
RETURN 
END 

SUBROUTINE  l'EUT(E  ,  T E # RHOE , DET , XL ) 

THIS  ROUTINE  CALCULATES  ENERGI  LOSS  TO  THE  T HER HAL  ELECTRCNS 
USING  FITS  PV  SCMwARTZ,  NISBET,  AND  GREEN  JGR  76,8425  (1971) 
-HERE  F:  IS  THE  PHOTOFLEC'IRON  ENERGX ,  TE  THE  THERMAL  ELECTRON 
TEMPERATURE,  RhOt  IS  THE  ELECTRON  DENS1TX,  CE/DT  IS  RETURNED 
VIA  DET,  L  VIA  XL. 

EE*«.618E-5*Tt 
X*(F-EE)/(fc-0.5J*EE) 

CETsU.O 
XL=0.0 

If  (X.LE.O.U)  GU  TO  l 
X«X**2,36 

DET*-2,E-44(RHCE**0,97)«X/(E**0,44) 

XI.*i,37E-I2*X/(  lRH0E**.03)*(E**U,94)  ) 

1  RETURN 
END 

0004 

(8( 1X.F9.2)  ) 
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-6 , lo6t-20  9,6221-19  1,0321-18  9.060E-19  4.273E-19  3.781E-19  1.706E-19 

1.302E-19  8.852E-20  4,2441-21  1,5821-21  3.55JE-22 

9.6221-19  1.0321-18  9,0601-19  8.5571-19  3.781E-19  3.403E-19  1.302E-19 

8.8521-20  9,9691-20  1,7651-21  3.5531-22  2.367E-22 

0 

3S5SU5)  10  9.15 

9.15  9.60  11.62  14.06  17,01  20.59  24,91  33,15  48,54  53.39 

-6.642E-19  2.959E-19  -8,5021-21  -2,0541-19  4,2691-20  4,2661-21  2.478E-21 

-1,1041-22  6.2531-22 

1.3281-18  -1.0271-20  -2.4841-19  5.1801-20  5,147fc-21  4.726E-21  -2,0611-22 

1.971E-22  -3.1271-22 

1.3451-19  -8.1151-19  4,1491-18  8.5671-18  3.364E-18  1.7031-18  3.293E-19 

1.5611-19  1,6221-20 

1.5091-18  4.992E-18  9,6751-18  4.295E-18  2.085E-18  8.609E-19  2.567E-19 

-3,6931-20  7.3551-21 

0 

0  4S3S3S  19  9.90 

9.50  12.00  14.50  17,00  19.50  22,00  24.50  29,50  37.00  44,50 

54,50  b2,0<>  74.50  82,00  102,00  127,00  152,00  172,00  199.50 

1,2321-20  1,2251-20  -3.03et-21  -3,3481-21  -3,6031-21  -2.846E-21  -1.U4E-21 

-4.1771-22  -1.6081-22  -6.1801-23  -2.6231-23  -6,7541-24  2.269E-24  1.431E-24 

1.7971-24  1.4481-24  1.396E-24  7,6721-25 

1 . 225E-20  -3.0361-21  -3.3481-21  -3.6031-21  -2.8461-21  -2.2281-21  -6,2651-22 

-1.808E-22  -8.2411-23  -1.9681-23  -1.1261-23  1,3611-24  3,8161-24  2,2461-24 
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1.448E-24 

1,1176-24 

1 , 055E-24 

8.185E-25 

•7 , 699E-20 

2,5846-19 

1 , 0S3E-18 

1 ,7346-18 

2.289E-18 

2,7091*18 

1 ,5326-18 

1.160E-18 

1.2306>«18 

9, 326E-19 

1,2126-18 

7.078E-19 

1,1 166-18 

4,0406-19 

2,9496-19 

2.660E-19 

3.035E-19 

2.062E-19 

2.5046-19 

1 , 053E-1 8 

1.734E-18 

2,2096-18 

2.709E-18 

3,0221-18 

1,7206-18 

1.230E-18 

1.240E-18 

9.097E-19 

1 , 178E-18 

6 , 695E- 1 9 

1 ,0796-18 

3,6916-19 

2.6606-19 

2.426E-19 

2.839E-19 

1 , 895E-19 

0  4S3P5P 

24  10.70 

10,70  13, 

20  15.70 

18,20  20, 

70  23,20  25 

,70  28.20 

30,70  33 

» 20 

35,70  38. 

20  40.70 

43,20  48, 

20  53,20  63 

.20  73.20 

85,70  100 

.70 

125.70  150. 

70  175,70 

198,20 

-1,6166-20 

-1 ,4336-20 

•1 ,5236-21 

-1 , 081E-21 

2 ,2446-23 

1 , 590E-2  2 

1.734E-22 

2,3556-22 

1.6476-22 

1,61 6E-22 

1 , 297E-22 

1 .069E-22 

9,5646-23 

3.689E-23 

2, 781E-23 

9.S3JE-24 

5, 347E-24 

2 , 542E-24 

1 » 192E-24 

3 , 866E-25 

1 ,6206-25 

8.669E-26 

5 . 057F-26 

-1.4336-20 

-1 .5236-21 

-1.081E-21 

2.244E-23 

1 ,5906-22 

1 ,7346-22 

2 , 355E-22 

1.6476-22 

1,6166-22 

1.297E-22 

1,0696-22 

9 , 564E-2  3 

7.379E-23 

2.781E-23 

1.907E-23 

5, 347F-24 

3 , 177E-24 

1.430E-24 

6.444E-25 

1.620E-25 

8,6696-26 

4.5516-26 

4,0126-26 

1.010E-19 

6.672E-19 

6,9596-19 

6.676E-19 

5,9876-19 

5,3066-19 

4.68SE-19 

4.129E-19 

3 . 662E- 1 9 

3.256E-19 

2.9UE-19 

2,6146-19 

2,3586-19 

1.062E-19 

8,8246-20 

3.654E-20 

2.709E-20 

1 , 653E-20 

1 .028E-20 

4,4306-21 

2,9756-21 

2.119E-21 

1 .770E-21 

b. 6726-19 

6,9596-19 

6.676E-19 

5, 987 E -19 

5, 306E-19 

4.6856-19 

4.129E-19 

3,6626-19 

3,2566-19 

2, 911 E- 19 

2.614E-19 

2 , 350E-  1  9 

2,1376-19 

8 . 824E-20 

7.450E-20 

2.7096-20 

2.084E-20 

1,2436-20 

7 , 655E-21 

2.975E-21 

2. 119E-21 

1 .5886-21 

1.4016-21 

C  4S3P3P 

22  11.00 

11.00  13. 

50  16,00 

18,50  21. 

00  23,50  26 

,00  28,50 

31,00  33 

,50 

36,00  38. 

50  43.50 

53.50  66, 

00  78.50  88 

,50  101,00 

126.00  151 

,00 

176.00  198. 

50 

-2.2V1E-20 

-4.488E-21 

-2 , 226E-22 

-4 , 459E-22  - 

8.661E-23  - 

5.088E-23 

-1.069E-23 

-7.703E-25 

1,3776-23 

9.692E-24 

1 ,7866-23 

7 .430E-24 

3 , 422E-24 

2.062E-24 

1. 3076-24 

1.1936-24 

6.781E-25 

2 , 599E-25 

1, 4516-25 

9.425E-26 

7,4 1 7E-26 

-4.4886-21 

-2 .2266-22 

-4.459E-22 

-8.661E-23  - 

5.888E-23  - 

1.0696-23 

•7 , 703E-25 

1 , 377E-23 

9,o92E-24 

1 .706E-23 

1.486E-23 

6.844E-24 

2.577E-24 

1.3076-24 

9.548E-25 

8.4766-25 

5. 199E-25 

1 , 45  IE-2  5 

9.425E-26 

6.6756-26 

5,0716-26 

1 ,4326-19 

3, 314F-19 

3.512E-19 

3.627E-19 

3.575E-19 

3.4916-19 

3.384E-19 

3,2746-19 

3,1636-19 

3.057E-19 

2 , 955E- 1 9 

1.428E-19 

6.681E-20 

4 , 760E-20 

4.217E-20 

4.753E-20 

3.527E-20 

1 .610E-20 

1 .401E-20 

1.2466-20 

1,2536-20 

3,3146-19 

3,5126-19 

3.627E-19 

3.575E-19 

3,4916-19 

3,3846-19 

3.274E-19 

3,1636-19 

3.057E-19 

2.955E-19 

2,8596-19 

1,3416-19 

5 , 964E-20 

4,2 1 7E-20 

3.797E-20 

4.414E-20 

3.245E-2C 

1.401E-20 

1 .246E-20 

1 .127E-20 

1 ,1576-20 

0  3030 

19  12.10 

12.10  14, 

6o  17.10 

19,60  22, 

10  24,60  27 

.10  29,60 

32.10  39 

,60 

47.10  54, 

60  64.60 

74,60  102, 

10  127.10  149 

,60  174,60 

199,60 
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1.00  /  1.00  «  HINTEREGGFR  FLUX.FITTFD  ATMOSPHF  RE  K t  0284U  .  7 

PHOTOE LFCTRP N  ENEPGY  GRID  IN  FV  -  5515.5 
100=NP  TFRRO 


•  00  00  OE 

00 

•  1000  OE 

01 

•  20  00  OF 

01 

•  3  00  0  OE 

01 

•  4000  OF 

0  1 

•5C000F 

01 

•  6000  OE 

01 

•  7 0 00  OE 

01 

•  6  0  00  OE 

01 

• 9COOOE 

01 

•  10  00  OE 

4  2 

•  1 10  00  E 

02 

•  1 2  on  OF 

02 

•  1 3000  E 

02 

. 14000E 

02 

•  1500  OE 

02 

.  1 6  00  OE 

02 

•  17  00  OF 

02 

•18000E 

02 

•  19  00  OE 

02 

•  2000  OF 

0  2 

•  21000E 

02 

.  22  00  OF 

02 

•2300CE 

02 

•  24  00  OE 

02 

.2 5  00  OE 

02 

•  2600  0  E 

02 

•  27  00  OE 

02 

•2800PE 

02 

•  29  00  OE 

02 

.3000  8F 

0  2 

. 21000E 

02 

.  32  00  OF 

02 

•33000E 

02 

•  34  00  OE 

02 

•  3500  OF 

02 

• 36000E 

02 

•  37  00  OE 

02 

•  ?  8000  E 

02 

•  3900  OE 

02 

•  4000  OE 

02 

•  41  CO  OE 

02 

•  4  2  0  0  0  E 

02 

•  4  3  000  E 

02 

• 44000F 

02 

•  4500  OE 

02 

•  46  00  OF 

02 

•47000E 

02 

•  4  8  00  C  E 

02 

• 45000E 

02 

•  5000  OF 

p  2 

• 51000E 

02 

•  52  00  OE 

02 

•5300CE 

02 

•  54  0  0  OE 

02 

.5500  OE 

0  2 

•  560  00  E 

02 

. 57000E 

02 

•  5  8000  E 

02 

.  5900  OE 

0? 

•  6000  OE 

p  2 

•  61 00  OF 

02 

•62000F 

02 

•  63  00  OE 

02 

•  64  00  OE 

0? 

•  65  0  0  OF 

02 

'•  6600  OE 

32 

•  670  00  F 

02 

•  68  00  0  F 

02 

• 69000E 

02 

•  7000  OE 

0  2 

•  71 00  OE 

02 

• 72000E 

02 

•  7  3  00  OE 

02 

•  74  00  0  E 

02 

.7500  OE 

o  2 

.•  76000E 

02 

• 77000F 

0? 

•  7  8  0  00  E 

02 

•  7500  OE 

02 

•  8000  OE 

o  2 

. 81000E 

02 

• 32000F 

02 

•83000E 

02 

.  8400  OF 

02 

•  8500  OF 

0  2 

•e6000E 

02 

•  87  00  OE 

02 

•8800CE 

02 

•  8900  OF 

02 

.9000  OF 

02 

•51300E 

02 

•  920  00  E 

C2 

•93000F 

02 

.  94000E 

02 

.9500  OF 

0  2 

• 96000E 

02 

o  97  00  OE 

02 

•98rOCE 

02 

•  9s  oo  or 

02 

SOLAR  ZENITH  53.649994 

ALTITUDE  4  09.000000 

TEMPERATURE  714.856934 

DENSITIES  HE  .523  7  0  OE  07 

N2  . 38 0  2 0  OF  06 

0  .21325  2C  08 

02  .241455F  04  25X.E16.6 

.ELECTRON  TcM  Pt  RA  TU'RE  .122000E  04 

ELECTRON  DENSITY  .241300E  06 

total  ion  nr  ns i t y  .ooooooe  oo 


PHOTOLLEJCTRPN  SPECTRUM 

6E12.4 

•  8  14 0 E  '00 

•6662 E  00 

. 45 17E  OC 

•1650E  OP 

• 2497E  00 

.  1  4  8  3  E  00 

•  292  IE  00 

.183*  E  OC 

. 1086E  00 

•1130E  00 

•3006E-01 

.428SE-01 

•1734E  00 

.3197  E-01 

•9165E- 02 

•  73  10  E-  01 

•2082E  00 

•5945E-01 

•2375E-01 

.211*1  E-01 

•9690E-01 

•  97  8  2E- 02 

•3993E  00 

•  53  03  E  OP 

.34435-01 

•2723E-31 

•5330E-01 

•4286E  00 

.71955-01 

•  68  7  8E -0  1 

•43115-01 

.9951 E-01 

• 4250E- 01 

•3647E-01 

•7718F-01 

.30  98E-01 

o  2046E -01 

.3069  E-01 

•1934E-01 

•1667E-01 

•2P93F-01 

.14  4ME-01 

.16215-01 

.2078  E-31 

•1959E-01 

•1778E-01 

•1607E-31 

•  31  ICE- P3 

•2164E-01 

.1697  E-01 

.  28  87E-0 1 

•2848E-01 

•2E41F-01 

•2936E-01 

•23245-01 

.3358  E-01 

.  1205E- 01 

•1673E-01 

•1699F-01 

•  2  6  32E -0? 

•7394E-02 

.2966  E-02 

•1234E-0? 

•1767E-0? 

•4726E-02 

•  69  7  6E  -  03 

•  3637  c  -'02 

.1444  r-02 

, 9774E-03 

•2444E-02 

•4343E-04 

•3562E-03 

•  3276 F -  03 

.2844  E- 04 

•699YE-05 

.1620E-03 

•7693E-04 

•2860E-03 

•  .21775-03 

.2343  E-34 

•2250E-03 

.  25  7 1  E-  03 

. 1059E-03 

.1757E-03 

•1317E-03 

•5568E-34 

. 8834E- 04 

. 17  9BE- 03 

.14 58 E-03 

.12 19E- 03 

•1549E-03 

•  1146  E- 33 

.  22  27  E- 03 

•9297E-04 

•1897E-03 

.15  60 E -  03 

•7661C-04 

.2068  E-03 

•9358E- 04 

•1153E-01 

SOLAR  ZENITH 

54,809953 

ALTITUDE 

384.330047 

TEMPFRATURE 

714.728271 

DENSITIES 

HE 

•  5569C  Of 

07 

N2 

•  102  JO  OF 

07 

0 

.358  1«7E 

04 

02 

• 7671C 9E 

04  25X.E16.6 

FHSC&Q1NG  PACK  BUNK  -  NO  T 

ELECTRON  TrMPERA  TL  RE 

.12  23  OOE 

04 

ELECTRON  CFNSITY  .306300T  36 
TuTAL  ICN  OFNFITY  .OOOOOOE  30 
PHOTOELFCTRPN  SPECTRUM  6E12.4 


•13825 

01 

•  1 1 16  E  01 

•8205E  00 

.28875 

OP 

.  4  3  09  E 

00 

•  25  7  0E 

on 

.49365 

00 

•  31 15  E  00 

•1853E  00 

•  17  57E 

00 

.5211E- 

01 

•7611c- 

n  1 

.P895F 

OC 

•54145-01 

• 1577C-0 l 

,  1  2  3  o  F 

p  n 

n  If  0  7  r 

0  P 

.  1  P  n  3  r 

0  n 

•  *40436-01 
•62666-01 

•  *71846-T)1 
•33965-01 
• 26806- 01 
•36436-01 
•39146-01 
•39966-02 
•61Z36-02 
•5451F-03 
#3875F-T)3 
•21716-03 
♦26046-03 
•12646-03 

SOUR  ZEN  I  TM 

ALTITUDE 

TETPFR*TU«E 

CfcNS I T  J  ES 


•33766-01 
•  •‘4777  E-01 

•  1678  £  00 

•  51 18  E-0 1 
•3473E-01 
•2852E-01 

•  5651 E -  01 
'  .4970E-Q2 

•  3109  £-02 

•  082  5E-04 
•42076-04 

•  .90406-04 
•1910 E-J3 

•  3065  6-03 


HE 

N2 

0 

02 


Electron  temperature 
Electron  density 

TtT>L  T  O'N  Dr  NETTY 
PHOTot  LFCTRp  N  SPECTRUM 

•  1007  E  01 

•  52506  oo 

•  9106  E -01 

•  5598  E  -0 1 
•8701 E-31 

•  281  fl E  00 

•  8098  E-01 

•  5787  E-01 
•0773E-01 

•  9028  E-01 

•  8  2 62  E  -Q  2 
•5191 E-02 

•  8139  E-00 

•  7687  E-00 
•1590E-03 
•3199E-03 
•5781 E-03 


• 22706  01 
'  «8293c  oo 

•  4799  c  '00 

•  7nilE-’01 

•  1197  c  00 
•l?00c  00 
.56505-01 
•00506-01 

•  a  1 0  2  6  -'0 1 

•65016-01 
■  .66216-02 
•1 02 06-01 
■  .90256-03 
•61526-03 
•3610E-03 
•0366E-03 
•21006-03 
SOLAR  ZENITH 
ALTITUDE 
TEMPFHATUifff 
DENSITIES 


HE 

N2 

0 

02 


ELECTRON  T6MiPE^ATU,RE 

electron  density 
1CT*l  ICIN  DENSITY 
PHOTOttf CTPPN  SPECTRUM 


•3889E  01 
•13506  01 
•76806 
•11916 
•22636  00 

•  1 °40  F  U  0 

•  9 1 15  f -o l 
•71526-01 
•9902E-01 
•1058E  ’00 

•  1 060  E  -  01 
•1600E 
-  •  «  k  u  r 


00 

00 


•01 

-  n? 


•  25  89  E  01 

•  35 8  0  E  00 

•  1085  E  00 

•  9027  E-01 

•  1581  E  00 

•05S6C  00 

•  1 360  E  00 
•9331 6-01 
•770CE-J1 

•  1522  E  00 

•  1328  E-31 
•9389 E-02 

•t'ldt. 1 < 


• 1589£  00 
•8796E-01 
•7152E-01 

•  31 79E- 01 
•3192E-01 
•0851E-D1 
•2021 E-01 
•2085E-02 
•1640E-02 
•12436-00 
• 3791E- 03 

•  10  5  3E- 03 
*3735 E-03 
•1507E-03 

55.990005 
360*69995 1 
710*096090 

•695S00E  07 
•  269 1 0  7E  07 
•590733E  08 
•233025C  OS 
•130000E  00 
•3809006  )6 

•  OOQOOOE  00 
6E1 2  »0 
•1358E  01 

•  31 03E  00 
•2795 E-01 
•2620E  00 
• 1 0  62E  00 
•150S£  00 
•52696-01 
•5280E-01 
• 907O£- 01 
•3373E-01 
•3O93E-02 
•2735E-02 
• 2238E-00 
•8331E-03 
•2420E-03 
•6226E-03 
• 2480E-03 

57.229  99  6 
338. 100098 
710.085893 

•  7 1 3  2  0  Qfc  07 

•  6  075  0  Of  07 
•950134E  03 
•683061 E  05 

•  1 3  65  7  OE  00 
•030900E  06 
•  0000006  00 
6El 2,4 
•21856  01 
•51  <J2E  OP 
•O878E-01 
•0175E  00 
•23536  00 

•  19  70  E  00 
•80506- 01 
•8061E-01 
•12996  00 
•50506-01 
•56706-02 

•  40  00  6-12 

•  tcnap-H 


•15916-01 
•71926  00 
•61 366-01 

.27336-01 

•29356-01 

•  47906- 01 
•2008E-01 
•29786-02 

•  91 00  E- 02 

•  26  90E-Q3 

•  42  61 E-03 
•30 18E-03 
•15066-03 
•1915E-0 1 


•6687E  00 
•12106  00 
•12116  00 
•3509E-01 
•25906-01 
•42816-01 
•28376-01 

•79506-32 

•61076-04 

•12556-03 

•17676-03 

•2439E-03 

•31426-03 


25X.E16.6 


•51 12E  00 
•28026  00 
•2056E  00 
.27136-01 
•11936  01 
•10266  00 
•45506-Oj 
•4870E-01 
•7972E-01 
•46 65 E-01 
•4975E-02 
•  64  076-02 
•44526-03 
•70 636-03 
•50396-03 
•2085E-03 
•3182E-01 


00 


•7061E 
•9101 F- 01 
•42996  00 

•  1 1 11 E  01 
•2026E  00 
•20176  00 
•58566-01 
•42876-01 

•  7 1 50E- 01 
•4699E-01 
•1325E-01 
•98106-04 
•2096E-03 
•29506-03 
•4P63E-03 
•52226-03 


25X,El6.ft 


■  »1 051E  01 

•  H56E  00 
•52126-01 

•  23  08E-01 
•509CE-01 
•0905c -01 

•  00  05E - 02 
•1180E-02 

•  60  02E-03 

•  48  02E-03 
»  29.626-03 
•2052E-0J 
•2528E-03 


«  99  00r  00 
.40396  OP 
•33196  00 
•45386-01 
•19106  01 
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CALCULATE  S  PRODUCTION  SPtCTRUM  OF  SECONDARY  ELECTRONS 
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DIMENSION  NORBI T(5) »PSS(200) 

TATA  NAMES/16H  HE  N  7  01  0?/ 

DATA  AFP/1. 443F-18.  1.134E-17.  4.750E-18.  9.692E-18/ 

CATA  S  IG/3 .48  OE -17  >  2.160E-16.  1.240E-16.  2.530E-16/ 

CATA  RHO/ 5 .37  OE  +  O  6 j  3.199F+08.  6.486F-.08.  1.698E+07/ 

DATA  WP/15.80  .  12.70  .  17.  40  .  17. <40  / 

PRINT  100 

100  FORMAT  (•  ENTER  PFLUX.N0R9IT  (  E  1  0  .3  .  5  X  .  5A  <4>  •  ) 

READ  (1.110)  PFLUX.NCRBIT 
110  FORMAT  (E10.3j5X.5A4) 

CO  120  K  =  1 .4 

P(K)*iRP0(K)»SI6(K)»PFLUX 
120  AR(K)=flHO(K)*AFP(K)»PeLUX 

DO  21 0  N= 1. 20 0 
ES=N— 0.5 
CC  =  0. 0 
CO  2C 0  K= 1.4 

2 on  CC=CC'+AR(K)/(  1.0*  (E  S/WP(* )  )  »»2.1) 

210  pcC(N)=CC 
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Abstract 

A  compilation  of  photoionization  and  photoabsorption 

cross  sections  is  presented  for  He,  0,  N2,  and  C>2  for  use  in 

studies  of  ion  and  photoelectron  production  in  the  terrestrial 

ionosphere.  In  wavelength  regions  where  rapid  variations 

occur  in  the  cross  sections,  averaged  cross  sections  are 

calculated.  When  necessary  the  cross  sections  have  been 

extrapolated  to  shorter  wavelengths.  The  cross  sections 

are  tabulated  at  the  wavelengths  of  the  solar  lines  and 

continua  given  in  the  solar  reference  spectrum  of  Hinteregger 
o  o 

from  M.030  A  to  ^34  A.  For  molecules,  N2  and  02 ,  branching  , 
ratios  are  given  for  ionization  into  the  ground  and  electronic 
states  of  the  molecular  ions  and  for  dissociative  ionization. 

* 

John  Simon  Guggenheim  Memorial  Foundation  Fellow  '78-' 79. 


I.  Introduction 


Any  detailed  theoretical  study  of  the  earth's  upper 
atmosphere  must  begin  with  a  calculation  of  the  production 
rates  for  major  ions  and  photoelectrons  due  to  solar 
radiation.  We  have  constructed  a  program  to  compute  these 

rates  using  a  model  atmosphere  or  measured  neutral  particle 

12  o 

densities  '  ,  the  solar  flux  tabulated  on  a  1  A  grid  and 

3  4 

observed  by  Hinteregger  ,  and  Heroux  and  Hinteregger  ,  and 

a  compilation  of  experimental  and  theoretical  photoionization 

and  photoabsorption  cross  sections  which  we  report  here. 

This  program  was  undertaken  in  conjunction  with  the  series 

of  Atmosphere  Explorer  satellites  (AE-C,D,E)  which  carried 

instruments  to  measure  simultaneously  many  geophysical 

parameters  of  the  terrestrial  ionosphere.  These  experiments 

have  been  described  in  volume  8  of  Radio  Science  (1973). 

The  neutral  species  of  interest  in  the  altitude  range  above 

120  km  are  He,  O,  N^,  and  0-,  with  photuionization  thresholds 
o  o  o  o 

at  504  A,  910  A,  796  A  and  1027  A,  respectively.  The  minor 
constituents,  N  and  NO,  are  not  included  in  this  compilation. 

We  have  been  guided  in  our  tabulation  by  the  intended 
atmospheric  applications  of  the  cross  sections,  for  which  a 
detailed  presentation  of  the  variaton  of  the  cross  sections 
is  neither  necessary  nor  desirable.  In  wavelength  regions 
in  which  the  cross  sections  are  highly  structured  we  have 
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frequently  computed  and  tabulated  averaged  values.  The 
accuracy  of  the  atmosphere  calculations  is  limited  here 
by  uncertainties  in  the  intensities  of  the  incident  solar 
radiation  and  in  the  densities  of  the  neutral  constituents. 
Many  of  the  measured  cross  sections  have  been  obtained  at 
low  spectral  resolution  and  at  room  temperature.  We  adopt 
these  nevertheless  but  the  possibility  of  important  errors 
should  be  noted’’ . 

Total  photoabsorption  cross  sections  for  the  atomic 
constituents  He  and  0  refer  only  to  photoionization.  For 
the  molecular  constituents  Nj  and  C>2  they  include  photoion¬ 
ization,  photodissociation  and  discrete  band  absorption  to 
excited  electronic  states.  In  photoionization  events, 
different  final  electronic  atomic  and  molecular  states  of 
the  ionic  products  may  be  populated  and  dissociative 
ionization  of  the  molecular  species  may  occur. 

All  the  photoionization  and  photoabsorption  cross 
sections  presented  here  have  been  tabulated,  using  interpola¬ 
tion  and  extrapolation  where  necessary,  at  the  solar  line 
and  continuum  wavelengths  of  the  Hinteregger  reference 
solar  flux^  which  is  based  on  measurements  by  the  Extreme 

Ultraviolet  Spectrophotometer  (EUVS)  on  board  AE-C^  and  on 
7  8 

rocket  data  '  .  We  restrict  the  tabulation  to  the  wavelength 
region  from  the  photoionization  threshold  of  07  at  1027  A 


to  33.74  A  which  is  the  shortest  wavelength  in  Hinteregger * s 

compilation  before  the  nitrogen  K-shell  ionization  edge  at 
o 

31  A.  The  intensity  of  the  solar  flux  decreases  rapidly 

o 

for  wavelengths  below  about  170  A  ,  so  that  the  solar  flux 

and  the  cross  section  data  at  short  wavelengths  are  not 

required  with  high  accuracy  for  most  aeronomy  applications. 

Specific  cross  sections  for  multiple  ionization,  are 

not  included  in  our  tabulation  except  that  the  total  photo- 

o 

absorption  cross  sections  include  them.  At  260  A  double 
ionization  of  contributes  only  about  2%  of  the  total 

9 

oscillator  strength  ,  and  although  the  multiple  ionization 

fraction  becomes  slightly  larger  at  shorter  wavelengths, 

the  solar  flux  is  decreasing  rapidly  in  this  region  and 

the  overall  effect  on  the  ion  abundance  is  negligible. 

Because  the  cross  section  data  have  been  obtained  from 

many  different  experimental  and  theoretical  sources  and 

have  involved  interpolation  and  extrapolation,  the  accuracy 

is  not  uniform.  Although  we  are  concerned  primarily  with 

valence  shell  rather  than  inner  shell  processes,  the  total 

cross  sections  of  the  molecular  species  at  short  wavelengths 
o 

(<150  A)  where  extrapolation  has  1'v.on  carried  out  include 
the  cross  sections  for  inner-shell  processes.  The  branching 

-f 

ratios  for  production  of  excited  electronic  states  of  ^ 
and  02+  are  unknown  in  the  wavelength  region  below  304  A, 
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and  lacking  any  data  we  have  used  the  304  A  values  down  to 
o 

34  A.  In  practice  the  branching  ratios  for  inner  shell 

absorption  presumably  increase  compared  to  the  valence 

shell  cross  sections  as  the  wavelength  decreases.  At 

o 

wavelengths  longer  than  ^650  A  there  are  several  regions 
in  the  molecular  cross  sections  which  are  densely  structured. 
Because  the  solar  flux  used  in  the  calculations  is  not 
tabulated  on  as  fine  a  wavelength  scale  as  the  oscillations 
in  the  cross  sections,  we  fitted  a  straight  line  through  the 
peaks  at  half  height.  In  other  regions,  in  which  the 
structure  was  less  dense  we  replaced  the  peaks  by  square 
waves  such  that  in  integrating  over  the  full  width  of  the 
lines,  using  che  grid  of  solar  wavelengths  given  by  Hinteregger 
the  effective  cross  section  is  equivalent  to  the  integrated 
cross  section  at  finer  resolution.  Thus  we  have  distorted 
the  shapes  and  magnitudes  of  the  cross  section  data  in  order 
to  obtain  equivalent  integrated  cross  sections  using  our 
solar  wavelength  scale.  At  particular  wavelengths  then, 
individual  cross  sections  may  be  incorrect.  In  addition,  there 
is  some  ambiguity  in  the  c^ross  sections  adopted  at  the  solar 
lines  due  to  the  limited  spacecraft  spectrometer  resolution, 
to  the  pressure  dependence  of  the  laboratory  data,  and  to  the 
uncertain  widths  of  the  solar  lines  compared  to  the  laboratory 
line  sources.  Thus,  errors  in  the  production  rates  could  occur 


if  there  exists  a  coincidence,  not  properly  included,  between 
a  resonance  for  an  absorbing  species  and  a  solar  line.  Where 
we  have  discovered  inconsistencies  between  several  sets  of 
data  we  have  generally  favored  the  recent  values.  The  following 
section  discusses  the  details  of  and  the  various  sources  for 
the  cross  sections  which  are  presented  in  Table  1.  Because 
cross  sections  for  each  species  are  treated  differently  in 
the  tables  due  to  the  diverse  nature  of  the  sources  for  the 
data,  it  is  desirable  to  consult  Section  II  before  using  any 
of  the  data  in  the  tables. 

We  have  incorporated  the  data  presented  here  and  the 

.  3 

tabulation  of  the  EUV  flux  given  by  Hinteregger  into  a  program 

which  calculates  the  ionization  rates  for  thermospheric 

constituents.  Table  A  shows  the  production  rate  in  s_1  for 

the  various  ion  states  in  an  atmosphere  with  zero  opacity. 

These  production  rates  are  valid  at  altiutudes  above  400 

kilometers  for  daytime  solar  zenith  angles. 

II .  Discussion  of  Cross  Section  Data 

o 

He:  For  He,  with  photoionization  threshold  at  504  A  , 

there  are  accurate  theoretical10'11  and  experimental12'13 

cross  sections  available  which  generally  agree  with  each 

other  to  within  about  5%.  The  cross  sections  near  threshold 
o  o 

(504  A  "  \  ^  454  A)  were  taken  from  Doyle,  Oppenheimer,  and 
Dalgarno14.  The  close-coupling  calculation  of  Jacobs10  for 
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a  simple  system  like  He  is  probably  accurate  to  about  2%, 

except  in  the  neighborhood  of  narrow  resonances.  We  converted 

Jacobs'  values1^  of  the  continuum  oscillator  strength  to  cross 

o 

sections  and  did  a  least  squares  fit  to  the  values  from  454  A 
o 

to  130  A.  As  the  cross  section  curve  is  entirely  smooth  in 

this  region,  little  error  should  result  from  this  procedure. 

o 

For  wavelengths  shorter  than  130  A  we  used  the  fit  to  extra¬ 
polate  the  cross  sections. 

O:  Photoionization  of  the  2p  valence  shell  electron  of 

atomic  oxygen  leads  to  0+(4S°),  0+(2D°),  and  0+(2P°)  with 

o  o  o 

thresholds  at  910.4  A,  732  A,  and  665  A,  respectively.  Vie 

used  the  empirical  fits  to  the  partial  cross  sections 

15 

for  these  three  channels  given  by  Henry  ,  normalized  to  the 

calculated  total  ionization  cross  sections  of  Taylor  and  Burke16 

Recent  measurements  near  threshold  by  'Kohl  et  al.17  tend  to 

confirm  the  values  of  Taylor  and  Burke  as  well  as  those  of 

1 8 

Pradhan  and  Saraph  ,  which  are  15-20%  higher  overall  than 

those  of  Henry1^.  Several  resonances  identified  by  Taylor 

and  Burke16  were  included,  from  706-608  A,  below  the  2D°  limit. 

Removing  the  2s  inner^shell  electron  of  atomic  oxygen 

gives  rise  to  0+(^pe)  and  0+(2Pe)  with  thresholds  at  435  A 
o 

and  315  A  respectively.  The  partial  cross  sections  for  this 
process  were  obtained  from  calculations  of  Dalgarno,  Henry 

19  2  0 

and  Stewart  as  modified  by  Henry  .  The  branching  ratios 

20  °  '21 

of  Henry  at  304  A  have  been  confirmed  by  Dchmer  and  Dehmer  . 
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Total  ionization  cross  sections  for  wavelengths  shorter  than 
o 

435  A  were  obtained  by  adding  these  partial  cross  sections 
for  inner  shell  ionization  to  the  total  ionization  cross 
sections  of  Taylor  and  Burke^. 

o 

N2 :  For  wavelengths  shorter  than  660  A,  the  ionization 

22 

efficiency  has  been  found  to  be  unity  ,  and  the  photoionization 

and  photoabsorption  cross  sections  are  equal.  In  the  region 
o  o 

180  A-650  A  we  interpolated  the  total  absorption  cross  section 

2  2 

data  of  Lee,  Carlson,  Judge  and  Ogawa  to  obtain  cross  sections 

at  each  wavelength  in  the  Hinteregger  reference  flux.  These 

authors  conservatively  estimated  the  error  in  the  cross  sections 

to  be  +20%,  and  as  there  is  almost  no  structure  in  this  region 

the  interpolated  cross  sections  should  be  of  the  same  accuracy. 

We  found  good  agreement  of  the  cross  sections  with  more  recent 

2  3 

data  of  Hamnet,  Stoll  and  Brion  as'  well  as  those  of  Gurtler, 

24 

Saile  and  Koch  .  We  have  not  included  the  peaks  attributed 

to  the  Rydberg  series  leading  to  the  state  of  N£+ 

°24 

between  500  and  550  A  and  they  are  probably  not  significant 

o  o 

for  our  studies.  From  180  A  to  34  A  the  data  were  extrapolated 

so  that  consistency  was  ^btained  with  the  absorption  cross 

oo.o  25 

sections  given  at  100  A,  68  A,  ard  44.6  A  by  Huffman  .  From 

o 

650  to  668  A  a  smooth  curve  was  drawn  joining  the  absorption 

22  c 

data  of  Lee  ct  al.  and  the  data  of  Huffman^,  passing  through 

<y  r 

two  points  measured  by  Samson,  Haddad  and  Gardner"  .  As  the 
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o 

ionization  decreases  from  100%  starting  at  660  A,  we  will 

now  discuss  absorption  and  ionization  cross  sections  separatel 

Where  considerable  structure  in  the  cross  section  is 
o  o 

evident,  from  668  A  to  734  A,  the  total  absorption  cross 

2  5 

sections  were  obtained  from  Huffman  .  We  estimated  a  back¬ 
ground  cross  section  by  drawing  in  a  base-line  on  the 
graphs  of  his  data;  superimposed  on  was  the  peak  cross 
section,  op.  Each  peak  area  was  approximated  as  a  square 
wave  over  the  same  wavelength  interval  as  the  actual  tri¬ 
angular  line-shape.  Any  dips  in  the  cross  sections  were 
approximated  in  similar  fashion.  The  total  cross  section 
at  each  wavelength  was  therefore  the  sum  (or  difference) 

O  O  o  n 

of  and  <jp.  From  734  A  to  986  A  Carter  has  tabulated 
oscillator  strengths  for  N2  between  adjacent  pairs  of  wave¬ 
lengths  and  again  a  square  shape  to  the  absorption  cross 

section  between  these  wavelengths  was  assumed.  Longward 
o 

of  986  A,  conflicting  measurements  for  the  N2  cross  section 

o 

have  been  reported,  and  above  1000  A  no  detectable  absorption 

2  8 

was  observed  by  Huffman,  Tanaka  and  Larrabee  .  Thus  from 
o  o 

986  A  to  1030  A  we  have  se^t  the  absorption  cross  sections 
to  zero. 

Photoionization  cross  sections  in  the  region  from 
660  ^  to  threshold  at  796  A  were  explicitly  calculated  using 
a  least  squares  fit  to  the  ionization  data  of  Cook  and 


109 


2  9 

Metzger  and  the  absorption  cross  sections  previously 

described.  Our  values  for  both  photoabsorption  and  photo- 

o 

ionization  cross  sections  in  the  wavelength  range  100  A 
o 

to  796  A  appear  to  be  reasonably  consistent  with  the  recent 

discrete  line  source  measurements  of  Samson,  Haddad  and 

26  9  o  o 

Gardner  and  the  work  of  Cole  and  Dexter  from  50  A  to  340  A. 

Photoionization  and  absorption  cross  sections  of  ^ 

25 

at  a  number  of  solar  lines  listed  in  Huffman's  Table  1 

were  included  explicitly.  The  majority  of  these  values 

30 

are  from  earlier  work  of  Samson  and  Cairns  .  Not  all  the 

25 

solar  lines  listed  in  Huffman  were  identical  with  those 

given  in  the  Hinteregger  reference  spectrum,  and  several 

lines  were  combined  and  the  cross  sections  averaged. 

The  dissociative  ionization  of  N^,  producing  N+  ions, 

was  treated  as  arising  from  a  single  state  with  threshold  at 
o 

509  A.  The  partial  cross  sections  for  this  process  were 
obtained  by  multiplying  the  total  ionization  cross  sections 
by  the  fractional  yield  for  dissociative  ionization  obtained 
from  Table  B.  Values  for  this  yield  were  derived  from  the 
data  of  Wight,  Van  der  Wiel  and  Brion^1  and  of  Fryar  and 
Browning  .  The  dissociative  ionization  yield,  Y,  can  be 
obtained  from  Table  B  by  linear  interpolation  between  the 

o 

listed  values,  with  the  exception  of  the  region  from  387  A 
to  477  A  where  the  following  form  should  be: 


1 10 


Y=0. 0329+8. 13xl0-6x(A-4 42) 2 

The  remaining  part  of  the  total  ionization  cross  section, 
c^on(l-Y),  is  apportioned  among  five  electronic  states  of  N2  + 
according  to  the  branching  ratios  listed  in  Table  C.  At  each 
wavelength,  the  branching  ratios  sum  to  1  and  at  wavelengths 
not  listed  in  Table  C  the  ratios  can  be  obtained  by  linear 
interpolation  between  adjacent  values.  Partial  photoionization 
cross  sections  are  obtained  by  multiplying  the  total  cross 
section  (minus  the  dissociative  ionization  cross  section) 
by  the  relevant  branching  ratio.  Branching  ratios  for  N2 

TO  OO  04 

photoionization  have  been  given  by  a  number  of  workers  '  '  ' 

and  the  agreement  among  the  published  results  is  generally 

within  10%.  There  appears  to  be  some  structure  present  in 

o  o 

the  published  branching  ratios  between  670  A  and  720  A,  but 

this  we  have  not  included  because  the  measurements  differ 

the  most  here  and  no  strong  solar  lines  appear  in  this  region. 

Since  no  measurements  have  been  reported  for  wavelengths 

o 

shorter  than  210  A,  we  assumed  constant  branching  ratios  from 
o  o  • 

210  A  to  34  A.  For  most  aeronomic. calculations  this  assumption 
is  of  minor  importance. 

02:  The  photoabsorption  and  photoionization  cross  sections 

4  o 

for  02  in  wavelength  range  612-34  A  were  obtained  in  an  analogous 
fashion  to  the  N2  cross  sections  in  this  renlon,  from  the  data 
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1 


of  Lee,  Carlson,  Judge  and  Ogawa^  and  from  the  short 

2  5 

wavelength  values  of  Huffman  .  The  ionization  yield  is 

o 

equal  to  unity  for  wavelengths  shorter  than  ^670  A.  For 

o 

wavelengths  longer  than  670  A,  a  great  deal  of  structure 
is  present  in  both  the  ionization  and  absorption  cross 
sections . 

o  o 

From  612  A  to  742  A  absorption  cross  sections  were 

obtained  by  graphically  interpolating  the  points  of  Samson, 

Gardner  and  Haddad^~*.  Such  a  fit  passes  roughly  through 

29 

the  highly  structured  region  of  Cook  and  Metzger  and  is 

in  general  agreement  with  the  overall  contour  of  their  data, 
o  o 

From  742  A  to  870  A  the  cross  sections  are  densely  structured 

and  we  used  a  least  squares  fit  to  the  absorption  data,  region 

by  region,  to  give  a  smooth  overall  profile  which  joins 
35  ° 

Samson's  data  at  742  A.  The  absorption  profile  in  the  870- 
o 

1030  A  region  is  marked  by  wider,  well-separated  peaks  and 

we  used  Huffman 's^^  data  to  fit  a  superposition  of  background 

and  peak  cross  sections,  and  a  as  described  for  N^.  Our 

o 

photoabsorption  cross  sections  for  from  50-350  A  are 

consistent  with  the  recent  measurements  of  Mehlman,  Ederer 

and  Saloman^,  Cole  and  Dexter^,  and  of  Samson  et  al.^"* 

o  o 

Photoionization  cross  sections  from  670  A  to  745  A 
were  obtained  by  a  smooth  fit  to  the  ionization  measurements 


of  Samson  et  al.  which  gives  a  curve  similar  in  shape  to 

\ 
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29 

the  overall  profile  of  the  earlier  Cook  and  Metzger  results. 

The  Samson  measurements  are  501-601  higher  than  the  Cook 

o 

and  Metzger  data  over  much  of  this  region;  from  74  5  A  lomiw.ud 

o  o 

the  two  sets  of  data  are  in  harmony.  From  745  A  to  870  A  we 

used  a  least  squares  fit  to  the  highly  structured  data  of 

o 

Cook  and  Metzger  and  in  the  region  870-1030  A  we  used  Cook 

•>9 

and  Metzger  ionization  data  and  total  absorption  cross 
sections  to  obtain  ionization  cross  sections  as  described 


for  N^.  At  the  solar  lines  listed  in  the  reference  sped. i urn 

we  adopted  the  ionization  and  absorption  cross  sections  listed 
25 

in  Huffman  ,  as  discussed  for  N2  . 

C' 

From  the  photionization  threshold  at  1027  A  \o  the 

o 

dissociative  ionization  threshold  at  662  A  we  have  used  the 

branching  ratios  for  production  of  C>2  +  electronic  states 

35 

given  by  Samson,  Gardner  and  Haddad  and  listed  in  Table  n. 

o  o  32 

In  the  region  304  A-662  A,  Fryar  and  Browning  have  measuied 
the  total  cross  section  for  dissociative  ionization  and 


obtained  values  which  exceed  the  sum  of  the  partial  cross 

sections  of  five  predissociating  02+  states  (labeled  4-8  in 

Table  D)  measured  by  Samson  et  al.'*'’  To  reconcile  this 

difference  we  include  an  additional  predissoeiat ing  branch 

o 

(labeled  0  in  Table  P)  with  threshold  at  (>(>2  A.  The  addition 
of  the  branch  made  renormalization  of  the  other  brandling 


ratios  necessary. 


As  no  measurements  have  been  reported  foi 


wavelengths  shorter  than  304  A  we  assumed  constant  branching 
o  o 

ratios  from  304  A  to  34  A.  Branching  ratios  at  any  wavelength 
o  o 

from  34  A  to  1027  A  can  be  obtained  by  linear  interpolation  of 
the  values  in  Table  D.  At  each  wavelength  branching  ratios 
sum  to  unity  and  the  partial  cross  sections  are  obtained  by 
multiplying  the  total  cross  section  by  the  relevant  branching 
ratio. 
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Explanation  of  Table  1. 


Photoionization  and  Photoabsorption  Cross  Sections 
for  0,  He,  N2  and  02 

- 1 8  2 

All  cross  sections  are  in  units  of  10  cm  . 

Before  using  the  table,  it  is  advisable  to  read  Section  II. 

LAMBDA  Wavelength  in  angstroms  for  which  the  solar  flux 

3 

is  given  by  Hinteregger  . 

0+(4S)  Cross  sections  for  photoionization  of  a  2p  valence 

0+ (2D)  electron  of  atomic  oxygen  giving  the  ground  state 

0+  ( 2 P )  0+(4S°)  or  excited  states  0+(2D°)  or  0+(2P°). 

0+ (4P)  Cross  sections  for  photoionization  of  a  2s  electron 

0+ ( 2 P * )  to  give  0+(4Pe)  or  0+(2Pe).  The  *  indicates  2Pe 

rather  than  2P°. 

T0T.0+  Oxygen  total  photoionization  cross  section. 

N2 ( ABS )  Total  photoabsorption  cross  sections  for  N2  and  02 

02(ABS)  respectively.  These  are  cross  sections  for 

ionization,  dissociation,  and  transitions  to 
excited  vibrational  and  rotational  levels. 

N2 ( I ON )  Total  photoionization  cross  sections  of  N2  and 
02 (10N)  respectively. 
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Table  A 


Ion  production  rates  in  s-^"  calculated  using  the  tabulated 

cross  sections  and  the  Hinteregger  solar  flux  from 
o  o 

1030  A  to  32  A  in  an  optically  thin  atmosphere 


Ion 

Rato 

He+ 

4. 84x10** 

oVs°> 

1.23x10*  7 

oVd°) 

9. 29xl0~8 

0+(2p°) 

5 . 42x10*" 8 

«2+<x2V’ 

3. 08x10*" 7 

N+  3 

3. 98xl0""8 

°2+«2V 

2.64xl0“7 

o2+(a4nu) 

1. 54xl0“7 

o+  b 

8.70xl0~8 

a)  N  produced 

by 

dissociative 

ionization 

of  Nr 

b)  0+  produced 

by 

dissociative 

ionization 

of  C>2  • 

Table  A:  In  constructing  the  table  it  has  been  assumed 

that  any  excited  atomic  and  bound  molecular  states  produced 

radiate  to  the  lowest  state  to  which  they  are  corrected  by  a 

dipole  transition.  For  example,  the  production  rates  for  all 

the  doublet  excited  states  of  are  included  in  the  production 

rate  of  0_+(X2II  ).  The  excited  0+(^PG)  state  decays  to 
i.  g 

while  the  0+(2Pe)  state  radiates  to  0+(~D°)  and  0+(2P°)  according 
to  the  branching  ratio  2.59:1^7.  The  ionization  caused  by 
absorption  of  these  "secondary  photons"  has  boon  nou looted. 
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Table  B 

Fractional  Yield  for  Dissociative  Ionization  of 


A  (A) 

Yield 

210 

0.360 

240 

0.346 

302 

0.202 

387* 

0.033 

477* 

0.041 

496 

0.024 

509 

0.000 

*See  text 

Table  B:  The  cross  section  for  dissociative  ionization 

,  00 

of  N2  at  any  waveleng^.i  from  32  A  to  509  A  is  obtained  by 

linear  interpolation  of  the  values  for  the  yield,  except  from 
o  o 

387  A  to  477  A  where  the  formula  in  Section  II  should  be  used, 
and  multiplying  by  the  total  ionization  cross  section. 
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Table  C 


Branching  Ratios  for  the  Photoionization  of 
See  text  for  sources 


Branch  1  2  3  45 

2  +■  2  2+  2  2  + 
D jsignation  X^E  A  n  aZ  PE  E 

g  u  u  u  g 


o 

X (A)  Branching  Ratios 


210 

0.271 

0.275 

0.110 

0.064 

0.278 

240 

0.271 

0.345 

0.110 

0.064 

0.210 

280 

0.271 

0.470 

0.095 

0.040 

0.124 

300 

0.271 

0.470 

0.110 

0.074 

0.075 

332 

0.300 

0.520 

0.120 

0.060 

0.000 

428 

0.460 

0.460 

0.080 

0.000 

500 

0.404 

0.506 

0.090 

600 

0.308 

0.589 

0.103 

660 

0.  308 

0.589 

0.103 

660.01 

0.308 

0.692 

0.000 

720 

0.420 

0.580 

747 

1.000 

0.000 

796 

1.000 

o  o 

Table  C:  From  796  A  to  509  A,  partial  photoionization  cross 

sections  can  be  obtained  by  multiplying  the  total  ionization  cross 

section  by  the  interpolated  branching  ratio.  For  wavelengths 

o 

shorter  than  509  A,  the  cross  section  due  to  dissociative  ionization 
must  first  be  subtracted  from  the  total  ionization  cross  section. 
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Branching  Ratios  for  the  Photoionization  of  02 
See  text  for  sources 


Branch 

1 

2 

3 

4 

5  6 

7 

8 

9 

Designation 

2 

XI 

g 

4  2 

an  +  a  n 

u  u 

b4Z  " 

g 

b2e  " 
g 

2-  - 
n  ci 

u  u 

2l  ' 
u 

2,4 1 

g 

o 

662  A 

o 

A  (A) 

Branching  Ratios 

304 

0.365 

0.205 

0.125 

0.055 

0.060  0.035 

0.030 

0.125 

0.000 

323 

0.374 

0.210 

0.124 

0.055 

0.060  0.035 

0.030 

0.000 

0.112 

454 

0.432 

0.243 

0.120 

0.055 

0.060  0.035 

0.000 

0.055 

461 

0.435 

0.245 

0.120 

0.055 

0.060  0.035 

0.050 

504 

0.384 

0.270 

0.126 

0.079 

0.026  0.000 

0.115 

537 

0.345 

0.290 

0.130 

0.098 

0.000 

0.137 

556 

0.356 

0.230 

0.225 

0.109 

0.080 

573 

0.365 

0.270 

0.216 

0.119 

0.030 

584 

0.306 

0.330 

0.210 

0.125 

0.030 

598 

0.230 

0.295 

0.375 

0.058 

0.045 

610 

0.235 

0.385 

0.305 

0.000 

0.075 

637 

0.245 

0.350 

0.370 

0.036 

645 

0.340 

0.305 

0.330 

0.025 

662 

0.270 

0.385 

0.345 

0.000 

634 

0.482 

0.518 

0.000 

704 

0.675 

0.325 

720 

0.565 

0.435 

737 

0.565 

0.435 

774 

1.000 

0.000 

1026 

1.00\) 

Table  D: 

For  molecular  oxygen 

,  the 

first  three  coliums  labeled  X":;  , 

a' 

a4':  +  A"  ,  ■ 

11  11 

4  - 

and  b  1  give  branching  ratios  to  bound  molecular  ion  states,  while  the  last  six 

columns  represent  dissociating  ionic  states.  The  sum  of  the  branching  rat ios 

o 

for  the  last  six  columns  for  wavelengths  less  than  662  A  mult  iplied  by  the  total  ion¬ 


ization  cross  section  «|ives  tin;  total  dissociat  i  ve  ioni.Mt  ion  cross  section. 
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